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Abstract 

In this paper, we analyze a PDE system arising in the modeling of phase transition and damage 
phenomena in thermoviscoelastic materials. The resulting evolution equations in the unknowns i? 
(absolute temperature), u (displacement), and X (phase/damage parameter) are strongly nonlinearly 
coupled. Moreover, the momentum equation for u contains X-dependent elliptic operators, which 
degenerate at the pure phases (corresponding to the values X = and X = 1), making the whole 
system degenerate. 

That is why, we have to resort to a suitable weak solvability notion for the analysis of the problem: 
it consists of the weak formulations of the heat and momentum equation, and, for the phase/damage 
parameter X, of a generalization of the principle of virtual powers, partially mutuated from the theory 
of rate-independent damage processes. 

To prove an existence result for this weak formulation, an approximating problem is introduced, 
where the elliptic degeneracy of the displacement equation is ruled out: in the framework of damage 
models, this corresponds to allowing for partial damage only. For such an approximate system, global- 
in-time existence and well-posedness results are established in various cases. Then, the passage to 
the limit to the degenerate system is performed via suitable variational techniques. 

Key words: Phase transitions, damage phenomena, thermoviscoelastic materials, elliptic degenerate 
operators, nonlocal operators, global existence of weak solutions, continuous dependence. 

AMS (MOS) subject classification: 35K65, 35K92, 35R11, 80A17, 74A45. 

1 Introduction 

We consider the following PDE system 

c(??)i?t+X t i? + p??div(ut)-div(K(^)V7?)=ff infix(0,T), (1.1) 
u tt - div(a(X)R u £(u t ) + b{X)R e e(u) - p&l) = f in fi x (0, T), (1.2) 

X t + ^(-oo,o] (X«) - div(d(x, VX)) + W'(X) 3 -b'(X) £ ( u )^ £ ( u ) + infix (0,T), (1.3) 

which describes a thermoviscoelastic system occupying a reference domain fi C R d , d € {2,3}, supple- 
mented with suitable initial and boundary conditions. The symbols i? and u respectively denote the 
absolute temperature of the system and the vector of small displacements. Depending on the choices of 
the functions a and b. we obtain a model 
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- for phase transitions: in this case, X is the order parameter, standing for the local proportion of 
one of the two phases; 

- for damage: in this case, X is the damage parameter, assessing the soundness of the material. 
We will assume that X takes values between and 1, choosing and 1 as reference values: 

- for the pure phases in phase change models (for example, X = stands for the solid phase and 
X = 1 for the liquid one in solid-liquid phase transitions, and one has < X < 1 in the so-called 
mushy regions); 

- for the completely damaged X = and the undamaged state X = 1, respectively, in damage models, 
while < X < 1 corresponds to partial damage. 



1.1 The model 

Let us now briefly illustrate the derivation of the PDE system (|l.l|) - (|1.3|l . We shall systematically refer 
for more details to [45], where we dealt with the case of phase transitions in thcrmoviscoelastic materials, 
and just underline here the main differences with respect to the discussion in |45j . 

Equation (|1.2[) , governing the evolution of the displacement u, is the classical balance equation for 
macroscopic movements (also known as the stress-strain relation), in which inertial effects are taken into 
account as well. It is derived from the principle of virtual power (cf. |18p. which yields 

u tt - div a = f in O x (0, T), (1.4) 

where the symbol div stands both for the scalar and for the vectorial divergence operator, a is the stress 
tensor, and f an exterior volume force. For a, we adopt the well-known constitutive law 

a = a nd + a d = —^ + —-^, 1.5 

with e(u) the linearized symmetric strain tensor, which in the (spatially) three-dimensional case is given 
by 6ij(u) := (iiij +Uj^)/2, i, j = 1,2,3 (with the commas we denote space derivatives). Hence, the 
explicit expression of a depends on the form of the free energy functional 5" and of the pseudopotential 
of dissipation 7. The former is a function of the state variables, namely X, its gradient VX, the absolute 
temperature and the linearized symmetric strain tensor e(u). According to Moreau's approach (cf. [T8] 
and references therein), we include dissipation in the model by means of the latter potential, which 
depends on the dissipative variables Vi?, X t , and e(u t ). We will make precise our choice for 5F and !P 
below, cf. (fl~T4")) and (fTTTj) . 

We shall supplement (|1.4[) with a zero Dirichlet boundary condition on the boundary of f2 

u = on c)fix(0,r), (1.6) 

yielding a pure displacement boundary value problem for u, according to the terminology of [12] . However, 
our analysis carries over to other kinds of boundary conditions on u, see Remark l2.9l 

Following Fremond's perspective, (|1.2|) is coupled with the equation of microscopic movements for 
the phase variable X (cf. [TS1 p. 5]), leading to (|1.3[) . Let B (a density of energy function) and H (an 
energy flux vector) represent the internal microscopic forces responsible for the mechanically induced heat 
sources, and let us denote by B d and H d their dissipative parts, and by B nd and H nd their non-dissipative 
parts. Standard constitutive relations yield 

B=B nd +B d = ^- + |^, (1.7) 

dX dX t 

H = H nd + H d = 41- + ■ (1-8) 

dVX dVX t y ' 
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Then, if the volume amount of mechanical energy provided to the domain by the external actions (which 
do not involve macroscopic motions) is zero, the equation for the microscopic motions can be written as 



S-divH = inOx(0,T), (1.9) 

where B and H will be specified according to the expression of 3 and CP. The natural boundary condition 
for this equation of motion is 

Hn = on<9ftx(0,T), 

where n is the outward unit normal to dfl. Thus (cf. (|1.14|Q we obtain the homogeneous Neumann 
boundary condition on X 

<9 n X = ondQx(0,T). (1.10) 

Finally, equation (jl.ip is derived from the internal energy balance 

e t + divq = g + a : e(u t ) + BX t + H ■ VX t infix(0,T), (1.11) 

where g denotes a heat source and e and q are obtained from 3 and CP by means of the standard 
constitutive relations 

<93 <9CP 

We couple equation (|l.lip with a no-flux boundary condition: 

q • n = on <90 x (0, T) 
implying (cf. (|1.17[l ) the homogeneous Neumann boundary condition 

<9 n tf = ondfix(0,T). (1.13) 



From the above relations and the following choices for the free energy functional and of the pseu- 
dopotcntial of dissipation (cf. (|1.14l) and (I1.17p ). we derive the PDE system (|1.1|) - (|1.3|) within the small 
perturbation assumption [23| (i.e. neglecting the quadratic terms \Xt\ 2 + a(X)e(u t )K v £(u t ) on the right- 
hand side of the heat equation). In agreement with Thermodynamics (cf. [T51[30] and [T^J Sec. 4, 6]), we 
choose the volumetric free energy 3 of the form 

3X0, e(u), X, VX) = J (f(&) + 6(x) £(u)I ^ eg(u) + 0( x , VX) + W(X) - - pt?tr(e(u))^ dx , (1.14) 

where / is a concave function of •&. Notice that the symmetric, positive-definite elasticity tensor R e is 
pre-multiplied by a function b of the phase/damage parameter X. In particular, 

- in the case of phase transitions in viscoelastic materials, a meaningful choice for b is b(X) = 1 — X, 
or a function vanishing at 1 (THl Sec. 4.5, pp. 42-43]. This reflects the fact that we have the full 
clastic contribution of 6(X)e(u)R e e(u) only in the non- viscous phase, and that such a contribution 
is null in the viscous one (i.e. when X = 1); 

- for damage models a significant choice is instead b(X) = X (cf. [50] and [TH1 Sec. 6.2, pp. 102-103] for 
further comments on this topic). The term ^ e ( u )^ £ ( u ) represents the classical elastic contribution 
in which the stiffness of the material decreases as X approaches 0, i.e. during the evolution of 
damage. 

The term <p(x, VX) + W(X) is a mixture or interaction free-energy. We shall suppose that <f> : f2 x M. d — > 
[0, +oo) is a normal integrand, such that for almost all x € VL the function (j>(x, •) : M. d — > [0, +oo) is convex, 
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C 1 , with p-growth, and p > d. Hence, the field d(x, ■) — W<f>(x, ■) : R d —> M. d , x € Q, leads to a p-Laplace 
type operator in ([L~5]) . The prototypical example is tf>(x, VX) = ||VX| P , yielding d(x, VX) := |VX| P " 2 VX. 
Let us point out that the gradient of X accounts for interfacial energy effects in phase transitions, and 
for the influence of damage at a material point, undamaged in its neighborhood, in damage models. In 
this sense we can say that the term i|VX| p models nonlocality of the phase transition or the damage 
process, i.e. the feature that a particular point is influenced by its surrounding. In damage, this leads 
to possible hardening or softening effects (cf. also [5] for further comments on this topic). Gradient 
regularizations of p-Laplacian type are often adopted in the mathematical papers on damage (see for 
example 1311033115), and in the modeling literature as well (cf., e.g., Q21 123 133 ) ■ In a different 

context, a p-Laplacian elliptic regularization with p > d has also been exploited in [T], in order to study a 
diffuse interface model for the flow of two viscous incompressible Newtonian fluids in a bounded domain. 

In the following, we will also scrutinize another kind of elliptic regularization in f| X . 3[) , given by the non- 
local s-Laplacian operator on the Sobolev-Slobodcckij space W S ' 2 (Q), hereafter denoted by A s (cf. (|2.35l) 
later on for its precise definition) . Recently, fractional Laplacian operators have been widely investigated 
(cf., e.g., [TTJ[52] and the references therein), and used in connection with real- world applications, such as 
thin obstacle problems, finance, material sciences, but also phase transition and damage phenomena (cf., 
e.g. [23 and [35)- For analytical reasons, we will have to assume s > d/2, which ensures the (compact) 
embedding W s - 2 (tt) <s C°(fi), in the same way as W 1,p (fl) <e C°(f2) for p > d. This property will play 
a crucial role in the degenerate limit to complete damage, as it did in [39| within the rate-independent 
context, cf. Remark 17.51 for more details. 

As for the potential W, we suppose that 

W = + 7, (1.15) 

with p : R — > (— oo, +oo] convex and possibly nonsmooth, and 7 : R — > R smooth and possibly nonconvex. 
We will take the domain of (3 to be contained in [0, 1]. Note that, in this way, the values outside [0, 1] 
(which indeed are not physically meaningful for the order parameter X, denoting a phase or damage 
proportion) are excluded. Typical examples of functionals which we can include in our analysis are the 
logarithmic potential 

W(r) =rln(r) + (l- r)ln(l-r)-cir 2 -c 2 r-C3 for r£ (0,1), (1.16) 

where c\ and Ci are positive constants, as well as the sum of the indicator function j3 := Imi] with a 
nonconvex 7. In such a case, in (|1.3[) the derivative W needs to be understood as the subdiffcrential 
dW = <9/3 + 7' in the sense of convex analysis. 

The term /M?tr(e(u)) in (| 1 . 14[) accounts for the thermal expansion of the system, with the thermal 
expansion coefficient p assumed to be constant (cf., e.g., |32p. Indeed, one could consider more general 
functions p depending, e.g., on the phase parameter X and vanishing when X = 0. This would be 
meaningful especially in damage models, where the terms associated with deformations should disappear 
once the material is completely damaged (cf., e.g., [5]). We will discuss the mathematical difficulties 
attached to this extension in Section 11.21 

For the pseudo-potential CP, following |T9j Sec. 4, 6] we take 

?(W, X t , e(ut)) - ^^|W| 2 + \\X t | 2 + m/^.o] (X t ) + a (x) £(Ut) y (Ut) , (1.17) 

where R„ is a symmetric and positive definite viscosity matrix, premultiplied by a function a of X. 
In particular, for phase change models, one can take for example a(X) = X. The underlying physical 
interpretation is that the viscosity term Xe(u t )K v e(u t ) vanishes when we are in the non-viscous phase, 
i.e. in the solid phase X = 0. Also in damage models the choice a(X) = X is considered, cf. e.g. [3T)] . 
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The heat conductivity function K will be assumed continuous; for the analysis of system (|1.1|) - (|1.3|) . 
we will need to impose some compatibility conditions on the growth of K(i?) and of the heat capacity 
function c(i3) = —#/"(#) in (jl.ip . see Hypothesis (II) in Section 12721 Furthermore, in (|1.17[) /i > is a 
non-ncgativc coefficient: for fi > we encompass in our model the unidirectionality constraint Xt < a.e. 
infix (0, T). In fact, throughout the paper we are going to use the term irreversible in connection with 
the case in which the process under consideration is unidirectional, which is indeed typical of damage 
phenomena. 

With straightforward computations, from (|1.7[t — (|1.9p and using the form of the free energy functional 
(|1.14p and of the pseudopotcntial of dissipation (|1.17[) . we derive equations (|l.l[) - (|1.3p . neglecting the 
quadratic contributions in the velocities on the right-hand side in (jl.ip by means of the aforementioned 
small perturbation assumption |23j . This is a simplification needed from the analytical point of view in 
order to solve the problem. Indeed, in a forthcoming paper we plan to tackle the PDE system (|l.ip - (|1.3p . 
featuring in addition these quadratic terms in the temperature equation. To do so, we are going to resort 
to specific techniques, partially mutuated from |17j . however confining the analysis to some particular 
cases. 

In fact, to our knowledge only few results are available on diffuse interface models in thermoviscoelas- 
ticity (i.e. also accounting for the evolution of the displacement variables, besides the temperature and 
the order parameter): among others, we quote [2TJ [22j E5] [46] . In all of these papers, the small per- 
turbation assumption is adopted. For, without it in the spatial three-dimensional case existence results 
seem to be out of reach, at the moment, even when the equation for displacements is neglected (whereas 
the existence of solutions to the full phase change model in the unknowns i? and X has been obtained in 
ID in [35]). This has led to the development of suitable weak solvability notions to handle (the usually 
neglected) quadratic terms, like in [T7] (where however u is still taken constant). Also in [15], a PDE 
system coupling the displacement and the temperature equation (with quadratic nonlinearities) and a 
rate-independent flow rule for an internal dissipative variable X (such as the damage parameter) has been 
analyzed. Rate-independence means that the evolution equation for X has no longer the gradient flow 
structure of f) 1 ,3[) : the term Xt therein is replaced by Sign(Xt), viz. in the pseudo-potential CP, instead 
of the quadratic contribution ^|Xt| 2 we have the 1-homogeneous dissipation term \X t \. In the frame of 
the (weak) energetic formulation for rate- independent systems |37j , suitably adapted to the temperature- 
dependent case, in [49| existence results have been obtained. A temperature-dependent, full model for 
(rate-dependent) damage has been addressed in [6] as well, with local-in-time existence results. 

1.2 Mathematical difficulties and related literature 

The main difficulties attached to the analysis of system (jl.ip - (|1.3p arc: 

1) the elliptic degeneracy of the momentum equation (|1.2p : in particular, we allow for the positive 
coefficients a(X) and b(X) to tend to zero simultaneously; 

2) the highly nonlinear coupling between the single equations, resulting in the the quadratic terms 
X t $, $div(u t ), and |e(u)| 2 in the heat and phase equations (jl.ip and (|1.3p . respectively; 

3) the poor regularity of the temperature variable, which brings about difficulties in dealing with the 
coupling between equations (jl.ip and (jl.2p when we consider the thermal expansion terms (i.e. we 
take p 0); 

4) the doubly nonlinear character of (|1.3p , due to the nonsmooth graph d/3 and the nonlinear operator 
— div(d(VX)) ~ — A p X (which on the other hand has a key regularizing role). Furthermore, if we 
set n > in (|1.3p to enforce an irreversible evolution for X, the simultaneous presence of the terms 
— div(d(VX)) and dI^ 00 fi](X t ) makes it difficult to derive suitable estimates for X, also due to the 
low regularity of the right-hand side of (|1.3p . 
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We now partially survey how each of these problems has been handled in the recent literature. 

As for 1), in [35JI35] we have focused on the phase transition case, in which a(X) = X and b(X) = 1 — X. 
We have proved the local-in time (in the 3-D-setting) and the global-in-time (in the ID-setting) well- 
posedness of a system in thermoviscoelasticity analogous to (|1.1[) - (|1.3[) (with the Laplacian instead of the 
p— Laplacian in (|1.3[) . in the case p = and p = 0, and for constant heat capacity and heat conductivity 
in (|l.l[l ). The main idea in [351135] to handle the possible elliptic degeneracy of p.2[) is in fact to prevent 
it. Specifically, we have shown that, if the initial datum Xq stays away from the values points and 
1, so does X during this evolution, guaranteeing that the operators in (|1.2[) are uniformly elliptic. This 
separation property is proved by exploiting a sufficient coercivity of W at the thresholds and 1, which 
for example holds true for the logarithmic potential 

In [7J [5] an isothermal (irreversible) model for damage has been considered: therein, because of the 
elliptic degeneracy of (|1.2p . the authors only prove a local-in-time existence result. For (isothermal) 
rate-independent damage models (38J, [H [39J [41] , the results change significantly: in this realm, only poor 
time-regularity of the solution component X is to be expected, because the 1-homogeneous dissipation 
contribution in X t to CP just ensures BV-estimates for the function t i— > X(x,t). That is why, one has 
to resort to the aforementioned notion of energetic solution [37] . in which no time-derivatives of X are 
featured. Therefore, this concept is very flexible for analysis, and has allowed for handling the (degenerate) 
case of complete damage in 9, 39] by means of a specially devised formulation we will refer to later. 

Concerning problem 2), as already mentioned existence results have been obtained in [17] for a full 
model of phase transitions (in the reversible case p = and for constant u), even featuring the term 
\X t \ 2 on the right-hand side of the temperature equation. Therein, a suitable notion of weak solution 
is addressed, consisting of the phase equation, coupled with a total energy balance and a weak entropy 
inequality, for which existence is proved by relying on an iterative regularization procedure. This technique 
cannot be applied to system l|l.l|) - (|1.3|) . Nonetheless, let us mention that a key assumption in [T7] is a 
suitable growth of the heat conductivity K. Following [HHHT], here we will combine it with conditions 
on the heat capacity coefficient c to handle the quadratic nonlinearities Xti) and $div(u t ) in (jl.ip . 

Due to the lack of "good" a priori estimates for t? mentioned in 3), we will not be able to encompass in 
our analysis the case of a non-constant thermal expansion coefficient p, e.g. p(X) = X, which would still be 
interesting for damage [5]. Indeed, such a choice would lead to an additional term of the type fiXt div(u) 
in the heat equation, which we would not be able to handle without resorting to further rcgularizations, 
and possibly proving only local-in-time existence results. Nonetheless, let us stress that, especially in 
case of phase transition phenomena, the choice of a constant p is quite reasonable (cf., e.g., |32j). 

As for 4), in [27] (dealing with Cahn-Hilliard systems coupled with elasticity and damage processes; 
see also [IE]), the authors have devised a weak formulation of fj 1 . 3[) (in the irreversible case p = 1) which 
has allowed them to circumvent its triply nonlinear character. Such a formulation strongly relies on the 
special choice /3(X) = I[o.+oc)(X) (which, joint with the irreversibility constraint, still ensures that X takes 
values in the meaningful interval [0, 1], provided that Xo G [0, 1]). It consists of a one-sided variational 
inequality (i.e. with test functions having a fixed sign), and of an energy inequality, see (|1.19|) later. 

1.3 Our results 

Unlike [ISJHS], here we shall not enforce separation of X from the threshold values and 1, and accordingly 
we will allow for general initial configurations of X. Then, it is not to be expected that either of the 
coefficients a(X) and b(X) stay away from 0, which results in the elliptic degeneracy of the displacement 
equation (|1.2|) . To handle it, we shall approximate system (]1.1[) - (|1.3[) with a non-degencrating one, where 
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we replace (|1.2j) with 

u tt -&v((a{X)+S)R v e(vit) + (b(X) + 6)R e e(u)-p&l)=f infix(0,T), for 5 > 0. (1.18) 

Let us note that, to rule out the degeneracy, it is sufficient to truncate away from zero only the coefficient 
a(X) of the viscous part of the elliptic operator in the momentum equation. However, for technical 
reasons which will become apparent in Section [7] (cf. Rmk. \77I\i , when addressing the asymptotic analysis 
as S 1 to the degenerate limit, we will need to truncate the coefficient b(X) as well, resulting in (|1.18[) . 
In the analysis of (|1.18|) . we will distinguish the cases p = and p ^ 0: let us stress that, in the latter, 
there is an additional coupling between the heat and the momentum balance equations, which needs to 
be carefully handled and indeed requires strengthening of some of our assumptions. Furthermore, to 
avoid overburdening the paper we will tackle the case p ^ only for the reversible system (i.e. with 
p = 0). More specifically, in Theorems [TJ @] and [5] we will establish global-in-time existence results for 
the non-degenerating system (jl.ll 11.181 II. 3j) with p = 0, both in the reversible and in the irreversible 
cases. We will work under quite general assumptions on c and K, basically requiring that c and K are 
bounded from below and above by the sum of a bounded function and function behaving like a small 
power of i? (cf. Hypotheses (I) and (II) in Sec. 12. 2[) . In Theorem[5J we will handle the case p ^ 0, p = 
and prove the existence of global solutions to p. 11 11.181 II. 3|) , under the more restrictive assumption that 
K is bounded from below and above by a function behaving like d 2+l/ , with < (d — 2)/(d + 2) < v < 1, 
cf. Hypothesis (VIII). A continuous dependence result, yielding uniqueness of solutions, for the non- 
degenerating isothermal reversible system, possibly with p ^ 0, will be given in Theorem [3] Finally, we 
will address the degenerate limit S 1 in Theorem [5] in a less general setting, in particular confining 
ourselves to the case p = 0. In what follows, we give more details on Thms. [HE] 

Our first main result Thm. Q] states the existence of solutions to system (jl.ll 11.181 II .3[) with p = 0, in 
the reversible case p = 0, with the heat equation (jl.lj) suitably reformulated by means of an enthalpy 
transformation (cf. Sec. 12. 2j) , switching from the temperature variable $ to the enthalpy w. Already 
in the proof of this global- in-time existence result, a key role is played by the aforementioned p- growth 
assumption on the function <j) (jl,14j) with p > d. In fact, it enables us to derive an estimate for X in 
I/°°(0, T; W Al ' p (f2)), which in turns allows for a suitable regularity estimate on the displacement variable 
u, leading to a global-in-time bound on the quadratic nonlinearity |e(u)| 2 on the right-hand side of (jl.3j) . 
For further details we refer to the proof of Thm. [TJ developed by passing to the limit in a carefully 
designed time-discretization scheme and exploiting BocCARDO&GALLOUET-type estimates on -d. 

Relying on the stronger Hyp. (VIII), in the case p ^ 0, p = we will obtain enhanced estimates on 
(the sequence, constructed by time discretization, approximating) ?9, cf. also Remark |2 . 1 01 later on. These 
bounds and the related enhanced convergences will enable us to handle the (passage to the limit in the 
time discretization of the) thermal expansion terms in (jl.ip and (|1.18p . In this way, we will conclude the 
proof of the existence Theorem [2] for system (jl.ll 11.181 1 1 . 3[) in the case p ^ 0. 

In the reversible and isothermal case, continuous dependence of the solutions on the initial and problem 
data is proved in Thm. [3] under a slightly more restrictive condition on the field (/>, which is however 
satisfied in the prototypical case of the p-Laplacian operator. 

As already mentioned, in the irreversible case p > a major difficulty in the analysis of system (jl.ll 
11.181 11.3(1 stems from the simultaneous presence in (jl.3[) of the multivalued operators <9/(-oo,o](X*) and 
/3(X) = d(3(X), (cf. (|1.15|0 . as well as of the p-Laplacian type operator — div(d(cc, VX)), which still has a 
key role in providing global- in-time estimates for |e(u)| 2 . To tackle this problem, following the approach 
of [57] we restrict to the yet meaningful case — I[o,+oo) an d consider a suitable weak formulation of 



7 



(|1.3|) . It consists (cf. Definition 12.131 later on) of the one-sided variational inequality 

J (x t (t)<p + d(.T, vx(*)) • + $(t)<p + y(x(t))v + &'(x(t)) £(u(t))I ^ g(u(t)) y> - 0(t) v ) dx > o 

for all v? G W 1,p (n) with 99 < 0, for a.a. t G (0, T), with X t < 0, £ G &T[o,+ao)(X), 

(1.19a) 

and of the following energy inequality for all t G (0, T], for s = 0, and for almost all < s < t: 
I |X t | 2 da:dr + f (<f>(x, VX(t))+W{X{t))) dx 

Jn Jn t (1.19b) 

< J (<f>(x,VX{s))+W(X(s)))dx + J X t ^-6 / (X) £(u) ^ i£(u) +#j dxdr. 

In Sec. 12.41 several comments and remarks shed light on this weak solvability notion for (| 1 . 3[) . In par- 
ticular, Proposition 12 . 141 shows that, if X is regular enough, (| 1 . and the subdifferential inclusion (|1.3I) 
are equivalent. In Theorem 0] we state the existence of global-in-time solutions to the weak formula- 
tion of system (jl.ll 11.181 II. 3[) with /j, > 0, consisting of the (weakly formulated) enthalpy equation, of 
(I1.18[) and of (|1.19|) . The proof is again carried out via a time-discrctization procedure, combined with 
Yosida-rcgularization techniques. 

Finally, Theorem [S] focuses on the isothermal case, i.e. with a fixed temperature profile. In this setting, 
we succeed in proving enhanced regularity for X, thus solving (|1.3[) in a stronger sense than (|1.19ll . In the 
particular case (f>{x, VX) = ^|VX| P the crucial estimate consists in testing (|1.3[) by dt(A p X + /3(X)) (where 
for simplicity we write /3 as single- valued) . This enables us to estimate separately the terms 9/(_oo 01 (Xt) 
(again written as single- valued) , A p X, and f3(X) in L°°(0, T; L 2 (£l)), which is the key step for proving the 
existence of solutions to the pointwise subdifferential inclusion (|1.3I) . 

Uniqueness results for the irreversible system, even in the isothermal case, do not seem to be at hand, 
due to the triply nonlinear character of equation 31) . cf. also Remark 12.181 ahead. Nonetheless, both in 
the reversible and in the irreversible case, in Thms. [TJ [5] and U wc will prove positivity of the temperature 
In fact, under suitable conditions on the initial temperature, for \i > wc will also obtain a strictly 
positive lower bound for 

For the analysis of the degenerate limit 5 i of p. 11 11.181 1 1 . 3[) . we have carefully adapted to the 
present setting techniques from [5] and |39j . These two papers deal with complete damage in the fully 
rate-independent case, and, respectively, for a system featuring a rate-independent damage flow rule for X 
and a displacement equation with viscosity and inertia according to Kelvin- Voigt rheology. In particular, 
we have extended the results from |39j to the case of a rate- dependent equation for X, also coupled with 
the temperature equation. Following the key observation is that, for any family (ws, U5, Xs)s of 

solutions to p. 11 11.181 II. 3[) (where w denotes the enthalpy), it is possible to deduce for the quantities 
fi s := y/a{X&) + Se(d t us) and rj s := y/b(Xs) + Se(ug) the estimates 

llA l ( 5ll-L 2 (0,T;L 2 (n;R d x d )): II Vs 1 1 L=° (0,T;L 2 (0;R d * d )) < C 

for a positive constant independent of S. Therefore, there exist \± and r\ such that, up to a subse- 
quence n s -± n in L 2 (0,T;L 2 (fl;R dxd )) and 77^*77 in L°°(0, T; L 2 (Q; R dxd )) as S I 0. According the 
terminology of [39], we refer to /_* and 77, respectively, as the viscous and elastic quasi- stresses. 

In Theorem [6] we will focus on the degenerate limit S I 0, confining the discussion to the case where 
p = and /j, > (viz. the map t i-> X(i, x) is nonincreasing for all x G fi). We refer to Remark 17.51 for 
a thorough justification of these choices. Passing to the limit as S J. in (| 1 . 18[) and exploiting the above 
convergences for (ng)s and (r] s )s we will prove that there exist a triple (u, /_t, rf) solving the generalized 
momentum balance 

u tt - div( V^POR^At + ^6(X)R e r7) = f in O x (0, T), (1.20a) 
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such that the quasi-strcsscs fulfill 



H = v / a(X)e(u t ), r) = \/b(X) e(u) a.e. in any open set A C n x (0,T) s.t. A C {X > 0}. (1.20b) 



In addition to (|1.20a[) - (|1.20b[) . the notion of weak solution to system (|1.1[) - ([1.3[) arising in the limit S I 
consists of the (weak formulation of the) enthalpy equation, of the one-sided variational inequality 



( (X t + 7(X)) <p + d{x, VX) • Vv») < 



1 



-jyRe t] + # I (/sdxdi 



v 26(X) ' - ' (1.20c) 
for all <p G £ p (0, T; W 1,p (fi)) n L°°(Q) with </? > and supp(p) C {X > 0}, 

and of a generalized total energy inequality, featuring the quasi-stresses rj and /x. While referring to 
Remark l7.4l for more comments in this direction, we may observe here that (|1.20c[) is in fact the integrated 
version in terms of quasi-stresses of the variational inequality (jl.19a.ll . 

Plan of the paper. In the next Section [2] we introduce the variational formulation for the initial 
boundary value problem associated to the PDE system (jl.ip - (jl.3[) . as well as our main assumptions. Then, 
we state Theorems [1] [5] on the existence/uniqueness of solutions for the reversible and the irreversible 
non- degenerating systems (i.e. S > 0). The existence Thms. [TJ [2J 01 and[5]rcly on the time-discretization 
procedure of Section [3] their proof is carried out by passing to the limit with the time discretization in 
Sections 14. 11 14. 21 16. 11 and 16.21 The continuous dependence Thm.[3]is proved in Section [5j Finally, Section 
[7] is devoted to the passage to the degenerate limit 8 \, 0. 
The following table summarizes our results 



Results 


p = 




p = 0, 6 > 


Theorem □ (Sec. [13]): 3 


Theorem H (Sec. [2U): 3 


p ^ 0, 6 > 


Theorem [21 (Sec. \2^i: 3 




p = 0, $ constant, 5 > 


Theorem 0] (Sec. 1231): 
uniqueness 


Theorem 151 (Sec. ED : im- 
proved regularity 


p 7^ 0, •& constant, S > 


Theorem 1 (Sec. 1231): 
uniqueness 


Theorem El (Sec. [23]): im- 
proved regularity 


p = 0,510 




Theorem! (Sec. E]) : 3 de- 
generate case 



2 Setup and results for the non-degenerating system 
2.1 Notation and preliminaries 

Notation 2.1. Throughout the paper, given a Banach space X we shall denote by || • \\x its norm, and 
use the symbol (•, -) x for the duality pairing between X' and X. 
Hereafter, wc shall suppose that 



n c 



d G {2,3} is a bounded connected domain, with C 2 -boundary dfl. 
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We will identify both L 2 (d) and L 2 (Q;M. d ) with their dual spaces, and denote by (-, •) the scalar product 
in R d , by (•, -) L 2 (n) both the scalar product in L 2 (Q), and in L 2 (fi; R d ), and by H&(Sl; R d ) and flg(fi; R d ) 
the spaces 

^(ft;^) :={vgff 1 (Sl;M <i ) : v = on dtt}, endowed with the norm ||v||^ 1(f2) := / e(v) : e(v) dx, 

Jn 

H$(n;R d ) := {v G H 2 (n;R d ) : v = on dCl}. 
For a, p > I we will use the notation 

W+' p (fl) := {C G W a ' p (fl) : ((x) > for a.a.x G fl} and analogously for W^' p (f2). (2.1) 
We standardly denote by 

A : H 1 ^) ff^fi)' the operator (/ht,u) H1(n) := ^ Vu • Vwdx 

and, for any w £ H (CI), by m(w) := (w, l)#i(Q) its mean value. 

Given a (separable) Banach space X, we will denote by BV([0, T];X) (by C^, eak ([0, T];X), respectively), 
the space of functions from [0, T] with values in X that are defined at every t € [0, T] and have bounded 
variation on [0,T] (and are weakly continuous on [0,T], rcsp.) 

Finally, throughout the paper we shall denote by the symbols c, c', C, C various positive constants 
depending only on known quantities and by vt (respectively vtt) or (whenever it turns out to be more 
convenient) dtv (respectively duv) the first (respectively) second partial derivatives with respect to time 
of a function v. 

Preliminaries of mathematical elasticity. In what follows, we shall assume the material to be 
homogeneous and isotropic, so that the elasticity matrix R e in equation (|1.3[) may be represented by 



R e e(u) = Aitr(e(u))l + 2A 2 e(u), 

where Ai,A 2 > are the so-called Lame constants and 1 is the identity matrix. In order to state 
the variational formulation of the initial-boundary value problem for (|1.1[) - (|1.3[) . we need to introduce 
the bilinear forms related to the X-dcpendcnt elliptic operators appearing in (|1.2[) . Hence, given a non- 
negative function r\ G L°°(Cl), let us consider the continuous bilinear symmetric forms a e \(r]-, •), a v i s (r)-, •) : 
H^(Cl;R d ) x H^(fl;R d ) -> K defined for all u,v G H^(Cl;R d ) by 

a c i(r/u,v) := (- div(r?R e £(u)), v) Hl , Q . Rd) = Ai / rj div(u) div(v) + 2A 2 V" / r/e ii (u)e 4i (v), 

Jn ■ „._i Jn 



i,j=l ' 



d p 

a vis (r/u,v) := (- div(^R t ,e(u)), v) ffl(0;Rd) = ^ / r)£ i:j e y (u)etf(v), 

■ ■ i J Si 



(2.2) 



where (iij) G R dxd is the viscosity matrix R v . Now, by Korn's inequality (see eg [TH Thm. 6.3-3]), 
the forms <x e i(?7-, •) and a v i s (r)-,-) are Hq (CI; K d )-elliptic and continuous. Namely, there exist constants 
Ci, C 2 > 0, only depending on Ai and A 2 , such that such that for all u, v G i7p(0;R d ) 

a c i(r/u,u) > mf (ry(x))Ci||u||^ 1(0) , a vis (r ? u,u) > inf (r)(x)) Ci||u||^i (n) , (2.3) 

|aei(?7u,v)| + |av is (r7U,v)| < C2|MU°°(n)||u||. ff i(n)||v|| ir i(n). (2.4) 

We shall denote by £(??•) : H&(Cl;R d ) -> H^^W 1 ) and V(ry) : H%(Cl;R d ) -> i?" 1 ^;!^) the linear 
operators associated with a c i(7]-, •) and a v i s (ji-, •), respectively, namely 

(£ (?yv) ,w) ffl(0 . Rd) := a e i(r/v,w), (V (r?v) , w) H1(fi;Rli) := a vis (?7v,w) for all v, w G Fo(ft;R d ). 

(2.5) 
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It can be checked via an approximation argument that the following regularity results hold: 

if 77 € L°°(d) and u e H^(Q-.R d ), then £ (r?u) , V (r?u) e iT^Q; K d ), (2.6a) 
if 776 W 1,d (n) and ue ff^(fi;R rf ), then £ (7711) , V (7711) e £ 2 (0; M d ). (2.6b) 

Remark 2.2 (The anisotropic inhomogeneous case). In fact, the calculations we will develop extend to 
the case of an anisotropic and inhomogeneous material, for which the elasticity and viscosity matrices R e 
and arc of the form R e = (gijkh) and R„ = (£ijkh), with functions 

9ijkh, km eC^fi), i,j,k, h =1,2,3, (2.7) 

satisfying the classical symmetry and ellipticity conditions (with the usual summation convention) 

Qijhh — Qjihh — Qkhiji ^ijkh — ^jihh — &khij , k^h — 1, 2, 3 , ^ 

3 Ci > : 9ijkhiijS,kh > Cl£ij£ij, £ijkh£,ij£,kh > Ci&j&j for all £y : = £jj , i, j = 1,2,3. 

Clearly, (|2.8[) ensures (|2.3[1 . whereas not only does (|2.7[) imply (|2.4p . but the (^-regularity also allows us 
to perform the third a priori estimate of Section 13.31 rigorously. 



In what follows we will use the following elliptic regularity result (see e.g. [TH Thm. 6. 3-. 6, p. 296], cf. 
also [H p. 260]): 

3C 3 ,C 4 >0 Vu6ff 2 (O;R d ): C 3 \\u\\ mm < || div(e(u))|| L2(n) < C 4 ||u|| ff2(n) . (2.9) 



Finally, in the weak formulation of the momentum equation (|1.2[) . besides V and £ we will also make 
use of the operator 

G p : L 2 {Q) H-\Q;M. d ) defined by (e P W^) H1 d) := -p f tfdiv(v)dx. (2.10) 

Useful inequalities. We recall the celebrated Gagliardo-Nirenberg inequality (cf. [321 P- 125]) in a 
particular case: for all r, q € [1, +00], and for all w € L 9 (f2) such that Vi> € 2/(0), there holds 

IMU-(n) < C G N|kl| H , 1 .. (f ,)lklll; ( %, with \ = (J - £) + (1 - ^ < < 1, (2.11) 

the positive constant Cqn depending only on c?, r, (7, 0. Combining the compact embedding 

, , . 00 if d = 2, 

H$(n;K. d ) <g W^< d - T '(0;R d ), with <T = ^ for all 77 > 0, (2.12) 

[6 if d = 3, 

(where for d = 2 we mean that fl£(n ; R d ) <g W /1,9 (fi;R d ) for all 1 < g < 00), with 03 Thm. 16.4, p. 
102], we have 

Vg>0 3C e >0 VueH%(n;R d ) : ||e(u)|| Lll .-„ (n) < g||u|| ff2(n) + C e ||u|| L2(n) . (2.13) 

We will also make use of the compact Sobolev embedding 

W hp (Ct) g C°(fi) for p > d, with d > 2. (2.14) 

We conclude with the following Poincare-type inequality (cf. Lemma 2.2]), with m(w) the mean value 
of w. 

Vq>0 3C q >0 yweH\n): \\\w\«w\\ m{a) <C q (\\V(\w\ q w)\\ L 2 {n) + \m(w)\' !+1 ). (2.15) 
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2.2 Assumptions and weak formulations 

We enlist below our basic assumptions on the functions c, K, W, d in system (|l.ip - (|1.3[) . 
Hypothesis (I). We suppose that 

the function c : [0, +00) — > [0, +00) is continuous, and 

2d ( 2 - 16 ) 

3<Ji>v> a > c > V0 G [0,+oo) : c (l+^) CT " 1 < c(0) < c 1 {l+"d) IJl - 1 . 

Hypothesis (II). We assume that 

the function K : [0, +00) — > (0, +00) is continuous and 
3c 2 , c 3 > W G [0, +00) : c 2 c(tf) < K(i?) < ca(c(t?) + 1) . 

Hypothesis (III). We require 

a G C^R), 6 G C 2 QR) arc such that a(x), b(x) > for all x G [0, 1]. (2.18) 
Hypothesis (IV). We suppose that the potential W in (|1.3[) is given by W = (3 + 7, where 



(2.17) 



dom(/3) = [0, 1] , /3 : dom(/3) — > M is proper, l.s.c, convex; (2-19) 
jeC 2 (R). (2.20) 

Hereafter, we shall denote by f3 = 9/3 the subdifferential of /3, and set 7 := 7'. 
Hypothesis (V). We require that there exists 



(2.21) 



a Caratheodory integrand : Q x M d — > [0, +00) such that for a. a. ir G f2 
the map </>(x, •) : R d -> [0, +00) is convex, with 0) = 0, and in C 1 (K d ) , 

and, setting d := : f2 x R d — > R d , the following coercivity and growth conditions hold true: 

3p>d, c 4 , C5 , C6 >0 fora.a.^0 VCeM d : j fjf 7 Sl-n ( 2 ' 22 ) 

[ |d(.T,C)| <c 6 (l + |C| p x ). 

A generalization of the p-Laplace operator. We now consider the realization in L 2 (il) of <fi, i.e. 

f Ld>(x,VX(x))dx if 0(-,VX(-)) G L 1 ^), 
$ : L 2 (0) [0,+oo], $(X):=^ Jnn ' K 11 Vi K)) y) (2.23) 

I +00 otherwise. 

Relying on [23 Thm. 2.5, p. 22], it is possible to prove that <J> is convex and lower semicontinuous on 
L 2 (fl), with domain £>($) := W hp (fl) (due to (|2T22]lh its subdifferential 9$ : L 2 (f7) =4 L 2 (ft) is a 
maximal monotone operator. In (|1.3[) we will take the elliptic operator 

B := <9$ : L 2 (0) =4 L 2 (ft). (2.24) 

Clearly, ® is a generalization of the p-Laplace operator. Note that 



dom(B) := {v G W 1,p (0) : sup 

weD($)\{0} 



d(x, Vv(x)) ■ Vw(x) dx 



IHU 2 (fi) < +00 > , 



and (cf. |H Ex. 2.4]) 



(Bu,u;)iS(n) = / d(x,Vu(x)) • Vw(x)da; for all w G £>($) . (2.25) 
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In order to obtain further regularity and uniqueness results for system (|1.1[) - (|1.3[) . we will have to 
assume that either of the following additional hypotheses holds true. 
Hypothesis (VI). We require that the function (j> fulfills the p-coercivity condition 

3c 7 >0 VienVC^eK 11 : (d(x,C) - d(x,rj),C ~ v) > c 7 \C-v\ p , ( 2 -26) 
and it is Lipschitz with respect to x, viz. 

3L>0 Vx,yeOVCeR d : \<j>{x, Q ~ <t>{y, 01 < L\x - y\{l + |CH- (2.27) 

Remark 2.3 (A regularity result). It was proved in [501 Thm. 2, Rmk. 3.5] that, if in addition to 
Hypothesis (V) the function <j> fulfills Hypothesis (VI), then 

dom(S) C W 1+<T ' p (fl) for all < a < -, and 

P (2.28) 
V0 <( r<- 3C a >0VveW 1+t7 - p (n) : \\v\\ w i + „, P{Q) < C a ||B(«)|| £ » C n). 

Hypothesis (VII). We assume that <f> complies with (|2.27p and with the following convexity requirement 
3c 8 >03k>0 Vs e Q VC,r7 £ R d , C ^ : D^(x, £)VV > c 8 (n + |C| 2 )I»7| 2 - (2.29) 

Remark 2.4. Assumptions (|2.27[) and (|2.29p guarantee the validity of the following inequality (cf. [3"0j 
for a proof) 

3c 9 >0 VxenVC,V£^ d ■ (d(x,C)-d(x, V ),C-v)>c 9 (n+\C\ + \v\) p - 2 \C~v\ 2 , (2.30) 
which will play a crucial role in the proof of Thm. [3J 
Example 2.5. The two p-Laplacian operators 

A p (X):=-div(\VX\ p - 1 VX), p > d, (2.31) 

A p (X) := -div((l + |VX| 2 ) p/2 ), p>d, (2.32) 

are clearly of the form (f2~2"5|) . and comply with ([2T2"2")| and ([2"1>^ - ([2~2T|) (cf. [H Ch. I, 4— (iii)]). 

Observe that (|2.29p is fulfilled by the p-Laplacian operator A p (|2.32[) . whereas for the degenerate 
operator A p (|2.31f) . inequality (|2.29[) holds with k = 0. 

A nonlocal alternative to the p-Laplacian operator. As done in [31], we could replace the p- 
Laplacian-type operator 23 (|2.25p in (|2.59p with a linear operator, with domain compactly embedded in 
C°(0). More precisely, as in [31] we could choose 

23 := A s : H s (il) — > H s (fl)* with s > ^. (2.33) 

In (|2.33[) , H s (il) denotes the Sobolev-Slobodeckij space W s ' 2 (tt), endowed with the inner product 

(zi,Z2)n>(n) ■= {z u z 2 )l^(q.) +a s (z 1 ,z 2 ), 

where 

a.(*,*) :=yj n ^^(.-i) d ^ ^ 2 - 34 ) 

Indeed, since d e {2, 3}, we may suppose that s € (1, 2). Then, we denote by A s : H s (rt) — > H S (Q)* the 
associated operator, viz. 

(A s X,w) Hs(n) := a s (X,w) for every X, w e iJ s (fi). (2.35) 
Observe that, for s > d/2 we have H s {n) <e C (T!). 
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Enthalpy transformation. We now reformulate PDE system (jl.ip - (jl.3|) in terms of the enthalpy w, 
related to the absolute temperature $ via 



w = h(ti) with h(r) := / c(s)ds. (2.36) 
Jo 

It follows from (|2 . 1 6[) that the function h is strictly increasing on [0, +oo). Thus, we are entitled to define 

Q(w) := < K(w) := . (2.37) 

10 iftu<0, c{Q{w)) 

In terms of the enthalpy w, the PDE system l|l.l|) - (|1.3|l rewrites as 

w t + X t e{w) + p<d(w) div(u t ) - div(K(w)Vw) = g in fi x (0, T), (2.38) 
u tt - div(a(X)R„e(u t ) + b(X)R e e(u) - p<d(w)l) = f in X (0, T), (2.39) 

X t + pdI { ^ Q] (X t )-div(d(x,\7X)) + W'(X)3 4'(X)^ + 9H infix (0,T), (2.40) 

supplemented with the initial and boundary conditions (where n denotes the outward unit normal to <9f2) 

w(0)=w , u(0) = uo, u t (0)=v o> X(0)=X o in SI, (2.41) 
d n w = 0, u = 0, d n X = on dVL x (0, T). (2.42) 

Remark 2.6. The enthalpy transformation (|2.36p was proposed in and further developed in [57], in 
order to deal with PDE systems where a quasilinear internal energy balance analogous to (|1.1[) is coupled 
with rate-independent processes. The advantage of this change of variables, is that the nonlinear term 
c(#)#t in (jl.ip is replaced by the linear contribution wt in (|2.38[) . We will exploit this fact, when proving 
the existence of solutions to (an approximation of) system (|2.38p ~ (|2.42p by means of a time-discretization 
scheme. 

For later use, let us observe that Hyp. (I) implies 

3 d , di>0 Vu;e[0,+oo): d l {w 1/ ' 71 - 1) < Q(w) < d (w 1/rT + 1) , (2.43) 
the map w <— > Q(w) is Lipschitz continuous. 

A straightforward consequence of the first of (|2.43l) is that for every s G (l,oo) 

3C S > g L x (0) : ||e(t«)|| L . (n) < C.(||«;||y.% cn) + 1). (2.44) 

Moreover, Hypotheses (I) and (II) entail 

3c>0 VweK: c 2 <iC(ro)<c. (2.45) 

Finally, in order to deal with the case p ^ 0, we will adopt the following further assumption, which 
we directly state in terms of the function K instead of K and c, in replacement of Hypothesis (II). 
Hypothesis (VIII). We require that the function K defined in (|2.37p (where c fulfills Hyp. (I) and 
K : [0, +oo) — > (0, +oo) is a continuous function) satisfies 

Bcio > 3q > Vw G [0, +oo) : K(w) = c w (w 2q + l) . (2.46) 

Indeed, we could slightly weaken (|2.46p by prescribing that K is bounded from below and above by two 
functions behaving like w 2q , and we have restricted to (|2.46p for simplicity only. Let us stress that, if 
(I2.46P holds, K is no longer bounded from above. 
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Problem and Cauchy data. We suppose that bulk force f and the heat source g fulfill 

f e L 2 (0,T;L 2 {n;R d )), (2.47) 
g G L^O. T; L 1 ^)) n L 2 (0, T; ff 1 ^)'). (2.48) 
and that the initial data comply with 

i?o G L ai (fl) whence w := h(& ) & L 1 ^) , (2.49) 
u a e H 2 (Q;R d ), v £ H^(fl;M. d ) , (2.50) 
XoGdom(3), p(X ) eL 1 ^). (2.51) 

Variational formulation of the non- degenerating system. We now consider the non-degenerate 
version of system (|2.38p - (|2.42p : as already mentioned in the introduction, to rule out the elliptic degener- 
acy of the momentum equation (j2.39[) . it is indeed sufficient to truncate away from zero only the coefficient 
a(X), cf. (|2.58[) below. In specifying the variational formulation of the initial-boundary value problem 
for the non-degenerate system, due to the O-homogeneity of the operator 3/(_oo.o] (|2.92j) we will just 
distinguish the two cases p = and p = 1. We mention in advance that the L r (Q, T; T4 7l,r (ri))-rcgularity 
for w derives from BocCARDO&GALLOUET-type estimates [5] on the enthalpy equation, combined with 
the Gagliardo-Nirenbcrg inequality (|2.11[) . We refer to the forthcoming Sec. 13.31 and to [35] f° r all details. 

Problem 2.7. Given S > 0, p G {0, 1}, find functions 

w G i r (0,T;^ 1 ' r (fi))nL oo (0,T;L 1 (fi)) nBV([0,T\;W iy (Q)*) for every 1 <r < 4^4, (2.52) 

a + 1 

u G ff^O, T; £T 2 (ft; R d )) n W 1,oo (0, T; H%(Q; R d )) n ff 2 (0, T; L 2 (rj; R d )), (2.53) 

X G L°°(0, T; W 1 *^)) n tf^O, T; L 2 (ft)), (2.54) 

fulfilling the initial conditions 

u(0, x) = uo(a;), u t (0, x) = vo(x) fora.a.xeQ, (2.55) 

X(0,x)=X (x) fora.a.xefl, (2.56) 

i/ie equations 



tp(i) w(t)(dx) — / / wiy9tdxds+ / / Xt@{w)ipdxds + p / / div(u t )Q(w)<pdxds 
Jo Jn Jo Jn Jo Jn 

/ K (w)S7 w\7 ip dx ds = / g<f + woip(0)dx 
o Jn Jo Jn Jn 



(2.57) 



for all Lp := C°([0,T]; W 1,r ' (0)) nW iy (0,T; L r ' '(O)) and /or olite (0,T], 
u tt + V((a(X) + 5)u t ) + £(6(X)u) + ep(e(tu)) = f mi? _1 (f2;M d ) a.e. in (0,T), (2.58) 

and the subdifferential inclusion 

X t + M^(-oo,o] (Xt) + S(X) + P(X) + 7 (X) 9 -b'(X) £ ( u ^ e£ ( u ) + e(w) inW^iny a.e. in (0,T). (2.59) 

Remark 2.8. Since to G BV([0, T]; W 1,r '(0)*), for all t G [0,T] one has io(i) G W iy (ft)* . Combining 
this with the fact that w G L°°(0, T; i 1 (il)), we have that w(t) is a Radon measure on for all t € [0, T], 
which justifies the notation in the first integral term on the left-hand side of (|2.57p . Moreover, let us 
note that the BV-regularity w.r.t. time of the absolute temperature, which is mainly due to the presence 
of quadratic nonlinearities X t O(w) and div(u t )0(w) in (|2.57[) . is quite natural for this kind of problems 
(cf., e.g. 0, PI], and @S]). 
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Remark 2.9. The proof of our results could be carried out with suitable modifications in the case of 
Neumann boundary conditions on u. as well. We would also be able to handle the case of Neumann 
conditions on a portion Tq of 9fi and Dirichlet conditions on T\ := <9fi \ To (\^o\, |Ti | > 0), provided 
that the closures of the sets Tq and Ti do not intersect. Indeed, without the latter geometric condition, 
the elliptic regularity results ensuring the (crucial) Hq (fi; M rf )-regularity of u may fail to hold, see [T2l 
Chap. VI, Sec. 6.3]. 

In what follows, we will refer to system (|2.57p - (|2.59p with p = (with p = 1, respectively), as the 
(non-degenerating) reversible full system (irreversible full system, resp). In both cases p = and p, = 1, 
we will call isothermal the (non-degenerating) system (|2.58[) - (|2.59p . where 0(u>) in (|2.59[) is replaced by 
a given temperature profile 6*. 

2.3 Global existence and uniqueness results for the reversible system 

Our first main result states the existence of a solution (w, u, X) to the reversible full system under 
Hypotheses (I)-(V); under the further Hypothesis (VI) we are able to obtain some enhanced regularity 
for X. Its proof will be developed in Section U by passing to the limit in the time-discretization scheme 
set up in Sec. [3] 

Theorem 1 (Global existence for the full system, p — 0, p = 0). Let p = 0, p = 0, and assume 
Hypotheses (I)-(V) and conditions (|2.47[) - (|2.5ip on the data f, g,$o, Uo, Vo, Xq,. Then, 
1. Problem\2. 7| admits a solution (w, u, X), such that there exists 



i G L 2 (0, T; L 2 {VL)) with £(x, t) G f3{X(x, t)) for a.a. (x, t) G x (0, T), fulfilling (2.60) 
X t + 3(X) + g + 7 (X) = -b\X) £ ( u )y (") + Q ( w ) a .e. in tt x (0,T). (2.61) 
Furthermore, (w,u, X) satisfies the total energy equality 



w{t){6x) + \ f \u t (t)\ 2 dx+ f f |X t | 2 dxdr 



2 



(2.62) 



a vis ((a(X) + S)u t , u t ) dr + ^-a cl (b(X(t))u(t) , u(t)) + *(X(t)) + / W(X(t)) dx 

= I w(s)(dx) + l [ \u t ( S )\ 2 dx + ^a cl (b(X( S ))u(s),u( S ))+<Z>(X( S ))+ f W(X(s))dx 
n 1 Jn 1 Jn 

j f • u t dxdr + j j g dxdr for all < s <t <T. 
s Jn Js Jn 

2. If, in addition, <j> complies with Hypothesis (VI), then there holds 

X G L 2 (0,T;W 1+rT ' p (Cl)) for all < a < -. (2.63) 

3. Suppose that 

g(x, t)>0 for a.a. (x, t) G Q X (0, T). (2.64) 
Then, w > a.e. in Q x (0,T), hence 

${x, t) := G(w(x, t)) > for a.a. (x, t) G fi x (0, T). (2.65) 

We are going to prove the energy equality (|2.62p by testing (|2.57p by y> = 1, (|2.58p by u t , (|2.6ip by X t , 
adding the resulting relations, integrating in time, and developing the calculations at the end of the proof 
of Thm.[T]in SecO 
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We now turn to the case when the thermal expansion coefficient p ^ 0. As previously mentioned, 
to prove existence of solutions we need to replace Hyp. (II) with Hyp. (VIII), which has a key role in 
deriving the enhanced regularity (|2.67[) for w, cf. Remark 12.101 later on. 

Theorem 2 (Global existence for the full system, // = 0, p ^ 0). Let p = 0, p ^ 0, and assume Hypotheses 
(I) and (III)-(V) and conditions (|2.47[) - (|2.51[) on the data f, g,iDo, u o, v 0i -^o>- Suppose moreover that 
Hypothesis (VIII) is satisfied (in place of Hypothesis (II)), and that 

w G L 2 (Sl). (2.66) 

Then, 

1. Problem \2.7\ admits a solution (w, u, X) fulfilling (|2.60p - (|2.61[) . such that w has the further regularity 

w G L 2 {0, T- H 1 ^)) n L 2(q+1) (0, T; L e{q+1) (SI)) n L°°(0, T; L 2 {Sl)) 

6o + 6 (2-67) 

nBV([0,T];W 2 - 8 (fi)') ^ s =^±|, 

and the weak formulation of the enthalpy equation holds in the form 

ip{t) w(t)(dx) — / / wipt dx ds + / / X t Q(w)p>dxds + p / / div(u f )Q(w)pdxds 
Jo Jn Jo Jn Jo Jn 

/ K(w)A(pdxds = / / gp + / WQip(0)dx 
'o Jn Jo Jn Jn 

for all test functions p G := C°([0, T]; W 2 ^{Sl)) n J3" x (0, T; i 6/5 (rj)) ; where X"(io) - 
is a primitive of K . Moreover, (w,u,X) complies with the total energy equality (|2.62j) . 

2. If, in addition, (f> complies with Hypothesis (VI), then the further regularity result (|2.63l) holds true. 

3. If in addition g complies with (|2.64[) , then (|2.65|) holds. 

Remark 2.10 (Outlook to the enhanced regularity (|2.67|l ). Let us justify the additional regularity (|2.67[) 
for w, by developing on a purely formal level, enhanced estimates on the enthalpy equation (|2.38p , based 
on the stronger Hypothesis (VIII). Indeed, we (formally) choose ip = w as a test function for (|2.57[) : 
re- integrating by parts in time and exploiting (|2.46[) we obtain for any t G (0,T): 



(2.68) 



w(t)\ 2 dx + c 10 [ f (\w\ 2q + l)\Vw\ 2 dxds 
n Jo Jn 

1 



2 



n Jo Jn Jo Jn 



<o / \w(fl)\ A dx+ / / \g\\w\dx+ / / (|X t | + H|div(u t )|)|eH|Hda;ds 



(2.69) 



Now, we observe that J Q \w\ 2q \\7w\ 2 dx = l/(q + l) 2 / n |V(|u>| 9 w)| 2 dx and that, due to the Poincare 
inequality ([2~T5|) and to the fact that w G L°°(0, T; L 1 (Sl)), there holds 

llkl^ll^(n)<C(HV(k|^)||^ ( n)+l). 

Therefore, taking into account the continuous embedding i^ 1 (SI) C L 6 (Sl), for the l.h.s. of (|2.69[) we have 
the lower bound 

f f (\w\ 2q + l) \Vw\ 2 dxds > c f (||V W || 2 L2(f2) + ds ~ C - ( 2 - 7 °) 

J J J 

Clearly, relying on (|2.48|) we can absorb the second term on the r.h.s. of (|2.69p into its left-hand side. 
On the other hand, using the fact that I := Xt + div(u t ) G L 2 (0,T; L 2 (Sl)) and taking into account the 
growth (|2.43p of O, we can estimate the last summand on the r.h.s. by 

do f f \£\{w 1+1 ^+l)dxds < C f f \£\{w q+1 +l)dxds 

Jo Jn Jo Jn j-r, 
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where the first inequality follows from the fact that 1 + 1/cr < (3d + 2)/ (2d) < q + 1 thanks to (|2.16l) 
and (|2.46|) . and g > is chosen sufficiently small, in such a way as to absorb J Q * IIHli^+i) (q) ds into the 
r.h.s. of (|2.70[) . Plugging (|2.70[) and (|2.71D into (|2.69[) we immediately deduce an estimate for w in the 
space w e L 2 (0,T;H 1 (n))nL 2iq+1 ^0,T;L 6 ^+ 1 \n))nL oo (0,T;L 2 (n)). Observe that, as a consequence 
of we L 2 (9 +1 )(0,T;L 6 (« +1 )(f2)), we have w 2q+1 e L 1+1 / e (0,T; L 3+3 / e (Q)), with e = 2q + 1. Therefore, 

A(£(iy)) is estimated in L 1+1 / e (0, T; W 2,S (Q)'), 

with s = (6g + 6)/(4g + 5) the conjugate exponent of 3 + 3/(2q + 1). A comparison in (|2.38[) entails an 
estimate for w t G L 1 ^^; W 2 ' 3 (Cl)'), and we conclude (|2"1)T)) . 

In fact, the BV-cstimate for u> (and accordingly, the regularity required of the test functions) could be 
slightly improved by resorting to refined interpolation arguments: however, to avoid overburdening this 
exposition we choose not to detail this point. 

To prove Thm. [21 we will need to combine the time-discretization procedure for system (|2.57p - (|2.59[) . 
with a truncation of the function K in the elliptic operator of (|2.57[) . cf. Problem 13.31 later on. Hence, 
in order to make the estimates developed in Rmk. 12.101 rigorous, we will have to pass to the limit in two 
phases, first with the time-step, and then with the truncation parameter, cf. the discussion in Sec. 14.21 

For the isothermal reversible system, in both cases p — and p ^ 0, we obtain a continuous dependence 
result, in particular yielding uniqueness of solutions, under the additional convexity property for </> in 
Hypothesis (VII). Indeed, the latter ensures the monotonicity inequality (|2.30|) for d, which is crucial for 
the continuous dependence estimate. We also need to restrict to the case in which a is constant. 

Theorem 3 (Continuous dependence on the data for the isothermal system, p = 0, p € K). Let p = 0. 
p 6 R. Assume that Hypotheses (III)-(V) and (VII) are satisfied, and, in addition, that 

the function a is constant. (2.72) 

Let (ii, Uq, Vq, Xq), i = 1,2, be two sets of data complying with (|2.47|) and (|2.50|) - (|2.5ip . and, accordingly, 
let (ui, Xi), i = 1,2, be the associated solutions on some [0, T] with fixed temperature profiles Q(wi) = Oj € 
L 2 (0, T; L 2 (tt)). Set M := max^i^ |ll u i!lff 1 (o.T;_f/g(a : K ti )) + l| • Then there exists a positive constant So, 
depending on M, 6, T, and such that 

||ui — U2|| H /i,^(0.T;L 2 (O;R d ))nH 1 (0,T;_f/ 1 (n:R d )) + 11^1 — X 2 \\ L°° (0,T;L 2 (H))nLP(0,T; W 1 -p (O)) 

< So(\\ul - ug|| H i (n .Rd) + ||vj - vl\\ L 2 {Q , m + \\X\ - XolU 2 (n) (2.73) 

+ ||fl — f2||i2(0,T;JJ- 1 (U)) + ||®1 — ©2||l 2 (0,T;L 2 (O))) • 

In particular, the isothermal reversible system admits a unique solution (u, X) . 

The proof is postponed to Section [5] 

Remark 2.11. Let us notice that, if we consider the s-Laplacian (|2.34[) instead of the p-Laplacian (|2.25|) . 
the continuous dependence result stated in Theorem [3] still holds true also without assumption (|2.72[) . In 
this case, for any two solutions X\ and Xi of the isothermal reversible system, (|2.73[) yields an estimate 
on — X2\\l 2 (o,t-,h s (n))- For further details, we refer to Remark |5. II at the end of Sec. [5] 

2.4 Global existence results for the irreversible system 

Heuristics for weak solutions. As mentioned in the introduction, the major problem in dealing with 
the subdifferential inclusion (|2.59[) in the case p — 1 is the simultaneous presence of the three nonlinear, 
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maximal monotone operators 9-f(_oo,ol! -^i an( i which need to be properly identified when passing to 
the limit in the time-discretization scheme we are going to set up in Section [3J We now discuss the 
attached difficulties on a formal level, treating f3 and 9-T(_oo,o] as single-valued. 

It would be possible to handle j3 = d(3 by exploiting the strong-weak closcdncss (in the sense of 
graphs) of the (induced operator) /3 : £ 2 (0, T; L 2 (tt)) =4 L 2 (0, T; L 2 (Cl)). Nonetheless, a L 2 (0, T; L 2 (n))- 
estimate of the term (3(X) in t|2.59|) cannot be obtained without estimating as well CB(X) (and hence 
dI(-oo,o}(Xt) by comparison), in L 2 (0, T; L 2 (f2)). To our knowledge, this can be proved by testing (|2.59D 
by d t (CB(X) + [3{X)) (cf. also [8]). The related calculations (which we will develop in Sec. [3l on the time- 
discrete level, for the isothermal irreversible system) would involve an integration by parts of the terms 
on the right-hand side of (|2.59|) . Thus, they would rely on an estimate in M /1 ' 1 (0, T; L 2 (Cl)) of the term 
— b'(X) e ( u )R<= £ ( u ) 4- Q(u;). However, presently this enhanced bound for &(w) does not seem to be at hand 
due to the poor time- regularity of w, cf. (|2.52D . 

That is why, for the temperature-dependent irreversible system we are only able to obtain the existence 
of solutions (w, u, X) to a suitable weak formulation of (|2.59p . mutuated from |27j . where we also restrict 
to the particular case in which 

P = I[0, +00 ). (2.74) 

Remark 2.12. In the present irreversible context it is sufficient to choose (3 as in (|2.74[) to enforce the 
constraint X £ [0, 1] a.e. in CI x (0,T). Indeed, starting from an initial datum Xo < 1 a.c. in CI we will 
obtain by irreversibility that X(-,t) < Xo < 1 a.e. in CI, for almost all t £ (0,T). 

The underlying motivation for the weak formulation of (|2.59[) we will consider is that, due to the 
1-homogeneity of I(-oo,o] , it is not difficult to check that (|2.59[) is equivalent to the system 

X t (x,t)<Q for a.a. (x,t) £ CI x (0,T), (2.75a) 

(X t (t) + 2(X(t)) + m + 7 (X(*)) + b'(X(t)) £ ( u (^))Re<u(t)) _ Q ^ Q 

2 r.pffi) (2.75b) 

for all <p £ W^ p (Cl) for a.a. t £ (0,T), 

(x t (t) + B(x(t)) + m + 7 (x(t)) + b'(x(t)) £ ( U (*))R* £ ( U (*)) _ Q (w(t)),X t (t)) <0 , 

2 r^(n) (2.75c) 

for a.a.t e (0,T), 

with £ £ <9/[ 0] +oo) (X) a.e. in 17 x (0, T) (and (•, ■)v^i,P(n) denoting the duality pairing between W 1,p (Cl)' 
and VF 1,P (J7), cf. Notation 12. II) . In order to see this, it is sufficient to subtract (|2.75c|) from (|2.75b[) . and 
use the definition of dI[ +oo y However, for reasons analogous to those mentioned in the above lines, 
the proof of (|2.75c|) is at the moment an open problem. Therefore, following |27j . in the forthcoming 
Definition 12.131 wc weakly formulate (|2.75[) by means of (|2.75a|l . (an integrated version of) (I2.75bp . and 
the energy inequality (|2.78[) below, in place of (|2.75cl) . 

Definition 2.13 (Weak solution to the (non-degenerating) irreversible full system). Let fj, = 1. We call 
a triple (w, u, X) as in (|2.52j) - (|2.54[) a weak solution to Problem \2.7\ if, besides fulfilling the weak enthalpy 
and momentum equations (|2.57[) - (|2.58[) . it satisfies X t (x,t) < for almost all (x,t) £ CI x (0, T), as well 
as 

(x t (t)<p + d(x, vx(t)) • V? + e(t)v + i(x(t))v + &'(x(t)) £(u(t)) ^ e£(u( ^ <p - e(w(t))v) dx > 

for all Lp £ Wb p (Cl), for a.a. t £ (0,T), 

(2.76) 

with £ £ 9/[o,+oo) (X) w following sense: 
CGL 1 (0 1 T;L 1 (n)) and (£(t), ^ - X(<)) r , (n) < e W^' p (tt), /or a.a. t G (0, T), (2.77) 
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and the energy inequality for all t G (0, T], for s = 0, and for almost all < s < t 

[ I |X t | 2 dxdr + $(X(i)) + [ W(X(t))dx 
J s Jn Jn 



t , , _r-, x ( 2 - 78 ) 



< $(X(s)) + / W{X{s))dx + I I X t I -b'{X) £ ^ R f^ + G{w) ) dxdr 



The following result sheds light on the properties of this solution concept. First of all, it states the total 
energy inequality (|2.T9[) : from (|2.79[) . we will deduce in Sec. [7] suitable estimates independent of 5, which 
will allow us to pass to the limit in Problem 12.71 as 5^0 for fx = 1. Furthermore, the second part of 
Proposition ^. 141 (whose proof closely follows the argument for Prop. 4.1]) shows that, if X is regular 
enough, then (|2.75ap and (|277S]> — ([2775]) are equivalent to (f2~59"l) . 

Proposition 2.14. Let (j, = 1. Then, any weak solution (w, u, X) in the sense of Dei. \2.13\ fulfills the 
total energy inequality for all t G (0,T], for s = 0, and for almost all < s < t 

w(t)(dx) + l f |u t (t)| 2 dx+ f f |X t | 2 dxdr 
^ Jn Js Jn 

a vis ((a(X) + S)ut, u t ) dr + ^-a cl {b{X{t))u{t),u{t j) + *(X(t)) + / W(X(t)) dx 

Jn (2.79) 
u>(*)(dx) + ^ / \u t (s)\ 2 dx+^a cl (b(X(s))u{s),u( S ))+^(X(s))+ I W(X{s))dx 

/ f ■ u t dxdr + / / g dxdr . 
/s Jo Js Jn' 

Assume now Hypotheses (III)-(V). Let (w,u,X) be as in (|2.52[) - ([2.54[) . and suppose in addition that 

B(X) G L 2 (0,T;L 2 (fi)) and t/iere exists £ G L 2 {0,T; L 2 (Q)) with 
£(x,t) G 5J [0i+oo) (X(a:,f)) /or a.a. (x,t) G SI x (0,T), 
swc/i i/iai (u;,u, X,£) comply with (|2.75ap and (|2~7d| - ((2~75|) . T/ien, 

3( G L 2 (0,T;L 2 (ft)) loiift C(M) € a/(_ oo , ](X t (x,i)) /or a.a. (x,t) G ft x (0,T) s.i. 
X t + ( + S(X) + £ + 7(X) = _6'(x)^l^M + ( w ) a . e , in fl x ( , T). 

Proof. In order to prove (|2.79p it is sufficient choose ip = 1 in (|2.57|l . test (|2.58p by u t , integrate in 
time, perform the calculations in the proof of Thm. [TJ and add the resulting equalities with the energy 
inequality (|2.78|) . The second part of the statement can be proved considering the energy functional 

S : L 2 {n) -> (-oo,+oo], <f(X):=$(X) + / W{X)dx. (2.82) 

Jn 

It follows from Hypotheses (IV) and (V), as well as the chain rule of [101 Lemma 3.3] that, if X complies 
with (|2.54p and (|2.80|) . then the map 1 1— > <§(X(t)) is absolutely continuous on (0,T) and fulfills 

d f 

—S(X(t))=j (%(X(x,t))+£(x,t)+j(X(x,t)))X t (x,t)dx for a.a. t G (0, T). (2.83) 
di Jo 

Therefore, differentiating (|2.78|l in time and using (|2.83p we conclude that (w, u, X, £) comply with (|2.75c|> . 
where the duality pairing (Xt,Xt) w i, P ^ is replaced by the scalar product in (Xt, X t )L 2 (n) ■ Likewise, 
(|2~75|) yields (|2.75b|) . Again on account of (j2~54) and (|2~%0"|) . it is not difficult to infer from (f2~75)) that 
_ b /( X ) £(u)R e e(u) + Q ^ _ ^ _ _ ^ _ ^ e a/ ( oo!0] (x t ) a .e. in ft x (0,T) as an element of 

i 2 (0,T;L 2 (n)), and (IMIl) follows. □ 
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Remark 2.15 (Energy inequality (|2.79|) vs. Energy identity (|2.62[) ). As we have already pointed out in 
the proof of Proposition ^. 141 to prove that any weak solution (w, u, X) in the sense of Dcf. 12.131 fulfills in 
the irreversible case (i.e. with p = 1) the total energy inequality (|2.79[) it is sufficient to choose ip = 1 in 
(|2.57p . test (|2.58l) by u t , integrate in time and add the resulting equalities to (|2.78|) . Instead, the proof 
of the total energy equality (|2.62p relies on the fact that, for p = we are able to obtain the stronger, 
pointwise form (|2.6ip of the subdiffcrential inclusion (|2.59[) . 

We can now state our existence result in the case p = 1, with p = 0. Its proof will be developed in Sec. 
16.11 by passing to the limit in the time-discretization scheme devised in Sec. [31 We mention in advance 
that, in this irreversible setting, in addition the basic assumptions of Hypotheses (I)-(V), we also have 
to require the p-coercivity condition (|2.26[) . It has a crucial role in proving strong convergence of the 
approximate solutions to X in L p (0, T; W 1,p (f2)), which enables us to obtain (|2.78p . Let us also highlight 
that, exploiting the additional feature of irreversibility, we prove a slightly more refined (in comparison 
with (|2.65[1 1 positivity result for the temperature cf. ()2.85p below. 

Theorem 4 (Existence of weak solutions for the full system, p = 1, p = 0). Let p = 1, p = 0, and 
assume Hypotheses (I)-(V) with (3 = /ro,+oo) as ^ n (|2 . 74|) . and conditions (|2.47p (|2.51l) on the data f, g, 
$o, Uo, vo, Xq. Suppose moreover that <fi complies with (j2.26p . Then, 

1. Problem \2.7\ admits a weak solution (w,u,X) (cf. Del \2.l3\) . 

2. Suppose in addition that 

g(x,t)>0 for a.a. (x,t) G O x (0,T) and 3i? > for a. a. x G : ■do(x)>tf o >0- (2.84) 
Then 

■&(x, t) := Q(w(x, t)) > ?? > for a.a. (x, t) G ft x (0, T). (2.85) 

Remark 2.16. If 23 is given by the s-Laplacian operator A s (|2.33p . the Definition 12. 131 of weak solution 
to the irreversible full system is obviously modified: in (|2.76p the term J n d(x, VX) ■ Wip is replaced by 
a s (X,tp) with tp G L 2 (0, T; Wl' 2 (Q)), whereas $(X) in (|2~75)l now reads \a s {X,X). Our existence result 
Theorem [4] extends to such a case, cf. also the forthcoming Remark l3.12l 

In Section [3 we will focus on the irreversible full system (|2.57p ~ (|2.59p with 23 = A s , and perform an 
asymptotic analysis of weak solutions (in the sense of Definition 12. 13[) in the degenerate limit 5 J, 0. 

Remark 2.17. Replacing Hyp. (II) with (VIII) and carefully tailoring the estimates and techniques for 
the proof of Thm. [5] (cf. Rmk. 12.1011 to the irreversible case, we could indeed prove the existence of weak 
(in the sense of Dcf. 12 . 1 3[) solutions also for p = 1 and p ^ 0. However we expect that, in the latter 
setting, only the weaker positivity result (|2.65p can be proved. Indeed, the estimate yielding the lower 
bound (|2.85p (cf. Step 4 in the proof of Lemma 13. 8p , cannot be performed on the enthalpy equation due 
to the presence of term — pQ(w) div(u t ). 

We finally turn to the irreversible isothermal case, and improve the existence result of Theorem 21 

Theorem 5 (Global existence for the isothermal system, p = 1). Let p = 1. In addition to Hypotheses 
(III)-(V), assume that 

b"(x)=0 for all x G [0, 1], (2.86) 

and suppose that the data f,Uo, Vq, Xq comply with conditions (|2.47p and (|2.50p - (|2.5ip . Suppose in 
addition fulfills (|2~2"B"| and (|2~27p . that 

23(X ), 0(X O ) e L 2 (n), (2.87) 
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and consider a fixed temperature profile 



6* e W 1 ' 1 (0,T;L 2 (a)). (2.88) 

Then, there exists a quadruple (u,X, £, £), fulfilling (|2.53[) ~ (|2.54[) . £ € /?(X) and £ € 9/(_ 00j o](Xt) a.e. 
m f2 x (0, T), as weH as 

X e L oo (0,r;W 1+CT ' p (O))n^ 1 ' oo (0,T;L 2 (fi)) for every < a < ^, (2.89) 

£e L°°(0,T;L 2 (n)), (2.90) 

Ce L°°(0,T;L 2 (O)), (2.91) 

satisfying equations (|2.58|) and (|2.81[) . ura£/i 0(ui) replaced by 0* . 



Remark 2.18. Uniqueness of solutions for the irreversible system, even in the isothermal case, is still 
an open problem. This is mainly due to the triply nonlinear character of (|2.59[) (cf. also [T3| for non- 
uniqueness examples for a general doubly nonlinear equation). 

A more general dissipation potential. As observed in Remark 16.11 later on, in Thm. [S] we could 
consider a more general dissipation potential in (|1.3[) . Indeed, in place of subdiffcrcntial operator dlt^^M , 
we could allow for a general cyclical monotone operator 

a := da : K =4 R, with 

(2.92) 

a : R — > M convex, lower semicontinuous, with dom(a) C (— oo, 0]. 

3 Time discretization 

First, in Section 13.11 we will approximate Problem 12.71 via time discretization. In fact, in the reversible 
case n = with p £ R, we will set up an implicit scheme (cf. Problems 13.21 and I3.3[) . whereas for the 
irreversible system fj, = 1 with p = 0, we will employ the semi-implicit scheme of Problem l3.4l Moreover, 
we will tackle separately the discretization of the isothermal irreversible system in Problem 13.61 We refer 
to Remarks 13.51 and 13.71 for a thorough comparison between the various time-discretization procedures, 
and more comments. Second, in Sec. 13.21 we will prove existence results for Problems I3.2ft3~6l Third, in 
Sec. I3.3l we will perform suitable a priori estimates. 

Notation 3.1. In what follows, also in view of the extension (|2.92[) mentioned at the end of Sec. 12.41 (cf. 
Rmk. IO]) . we will use a and a as place-holders for /(_oo.o] and dl^^m. 

3.1 Setup of the time discretization 

We consider an equidistant partition of [0, T], with time-step r > and nodes t^. := fcr, k = 0, . . . , K T . 
In this framework, we approximate the data f and g by local means, i.e. setting for all k — 1, . . . , K T 

f":=- t f(s)ds, g k T :=- t g(s)ds. (3.1) 

Problem 3.2 (Time discretization of the full reversible system, p, = 0, p = 0). Given 

w" := w , u" := u , u^ 1 := u - tv , X° := X , (3.2) 
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find {u^u*,**,^}^ C H 1 ^) x H${Q;R d ) x W^ity x L 2 {Q), with G /3(X^) a.e. in 0, fulfilling 

^6(^) + 4 *-i(u£) = fl * inif 1 ^)', (3.3) 



r r 



+ V ((a(X*) + <5) Ut " t ^ + £ (&(X?)u£) = f* a.c. in ft, (3.4) 

X " + 03(X*) + ^ + li^r) = _^ (x fc-i ) £ ( u r" 1 )y( u ^" 1 ) + e^) a . e . in 0, (3.5) 
r z 

where in (|3.3p the operator ^4 *-i : i? 1 (il) — > i7 1 (O)' is defined by 



(vl *-i(w),w) m . •= / Jjr(iw*- 1 )Vti; • Vvdx for all tu,« G H\n). 



(3.6) 



For the full reversible system with /? ^ 0, we work under the stronger Hypothesis (VIII) and thus 
prescribe a suitable growth on the function K, which is no longer bounded. Therefore, in order to 
properly deal with the elliptic operator in the enthalpy equation on the time-discrete level, we need to 
truncate K. We thus introduce the operator 

(A w *-i >M {w),v) HHQ) := J^K M (w k - l )Vw ■ Vvdx for all w,v G fl^fi) , (3.7) 



nth 



K{—M) if r < —M, 
K{r) if |r| < M, (3.8) 
K(M) if r > M. 



Observe that 



K M (r) > cio for every r G K (3.9) 
(with cio from (|2.46[) X Accordingly, we will have to truncate the function 0, replacing it with 

{0(-M) if r<-M, 
0(r) if |r| < M, (3.10) 
0(M) if r>M. 

Problem 3.3 (Time discretization of the full reversible system, fjb = 0, p ^ 0). Starting from (u°, u^T 1 , 
X°, io°) as QQ), /?nd u*, x£, C^}*^ G ff^fi) x ff£(ft;R rf ) x I^ 1 ^) x L 2 (ft) riA G /?(x£) 
a.e. in f2, fulfilling 



/c — 1 / k k—1 \ 

T —Q M (w k T ) + pdiv ( Ut ~" t J Q M {w k T ) + A w ,- l M (w k T ) = g k T m H 1 ^)', 



(3.11) 

u k - 2u k - 1 + ujT 2 „„ (, „ u$ - u^ -1 



k—1 / A; fc 1 \ 

^— + V f (a(X*) + S) " r j + £ (6(X?)u£) + e p {Q M {w L T )) = f r fc a.e. mfi, 

(3.12) 

+23(X^)+C fc +7(^) = -fc , (^ 1 ) £(U ^ 1)R 1 e£(U ^ 1) +6 M K fc ) a.e. mfi. (3.13) 
r 2 

We now present the time discretization of the full irreversible system, postponing to Remark 13.51 a 
detailed comparison between Problem 13.21 and the forthcoming Problem 13.41 Let us only mention in 
advance that, in the irreversible case we will restrict to the particular choice f3 = <9i]o.+oo)- Furthermore, 
in Problem [331 instead of the time discretization of (|2.59[) . we will consider the minimum problem (|3.19[) . 
such that its Euler equation is (|2.59[) discretized. We resort to this approach in view of the passage to 
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the limit argument as r — > 0, mutuated from [27| . which we will develop in the proof of Thm. Finally, 
due to technical reasons related to the proof of the Third a priori estimate in Sec. 13.31 we will also need 
to approximate the initial datum wq with a sequence 



(wot)t C W 1,r (£l) such that supr||Vwo T ||2f=(Q) < C, u>q t — > wq in asr^O 



r>0 



(3.14) 



with f = (d + 2)/(d + 1), cf. (|2.52[) . We construct (wq t ) t in such a way that, if #o complies with (|2.84|) . 
then for every r > 

W0r(a ; ) > w := /i(#o) > for a. a. x £ O. (3.15) 

Problem 3.4 (Time discretization of the irreversible full system, \i = 1, p = 0). Starting from the data 
(u°, u"\ X°, u;°) as in with w° = w 0t as in (l3~T4l . find {«;*, u*, X*, Cr }f=i G x #o (0| Rd ) x 

W Al -P(r2) x L 2 (r2), such that for all k = 1, . . . , # T there holds 



a.e. in 51 and 



a.e. in tt, 



and (<,u«,X*) fulfill 



x£ - X^ 1 



^eK^ 1 ) + A k-i (w k T ) = g k T in H 1 ^)', 



u£ - 2u 



fe-i 



,fc-2 



(a(Xj!)+<5)- 



£ (6(Xj)u* 



a.e. in SX 



(3.16) 

(3.17) 
(3.18) 



and 



X h T G Argmin^ 




X-Xl 



fc-i 



X - X 



fe-1 



dx + $(X)+ / (/3(X) +7(X))dx 



b/(x fc_ 1} £(u^)Re£(u^) _ Q^fc-l} 



(3.19) 

5} and ([3TT7]) - ([3TTg]l consists in the dis- 



Remark 3.5. The main difference between systems 
cretization of the coupling term 0(u>) in the temperature and in the phase parameter equations. Indeed, 
in the reversible case Q(w) on the l.h.s. of (|3.3p (and accordingly its coupled term on the r.h.s. of (|3.5|) . 
which will cancel out in the First a priori estimate of Sec. 13. 3[) . is kept implicit. Only relying on this 
we can prove the positivity of the discrete enthalpy w^, which for system (|3.3p - (|3.5[) would not follow 
from other considerations. Instead, in the time discretization (|3.17|1 - (|3.19|) we can allow for an explicit 
coupling term Q(w k ~ 1 ) in (|3.17|) and in (|3.19p . Therein, the positivity of the discrete enthalpy will be 
proved by means of a suitable test of the discrete enthalpy equation, relying on the irreversibility p. Kip . 

Because of its implicit character, in Lemma 13.91 existence for system (j3.3p ~ (|3.5p will be proved by 
resorting to fixed-point type existence results for elliptic systems featuring pseudo-monotone (cf. e.g. [3H1 
Chap. II]) operators. 

Instead, in the semi-implicit scheme of Problem l3.41 equations Q3.17p . (|3 . 18[) . and (|3 . 19[) are decoupled, 
hence we will proceed by tackling them separately, solving time-incremental minimization problems. Such 
a procedure could be useful for the numerical analysis of the problem. That is why, in Sec. 13.21 we will 
focus on the proof of Lemma l3~81 and develop in detail the calculations for system (|3.17p - (|3.19p . whereas 
we will only outline the argument for the existence of solutions to (|3.3|1 — (|3.5[) in Lemma 1331 



In the time-discretization of the irreversible isothermal system, we approximate the given temperature 
profile 6* (cf. (|2.88p ) by local means as well, i.e. 



6*" := - 



9*(s)ds. 



(3.20) 
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Contrary to the temperature-dependent irreversible case, we may again address a general maximal mono- 
tone P : R =4 K- However, in order to perform enhanced estimates on the discrete equation for X (cf. the 
Seventh and Eighth a priori estimates of Sec. 13. 3p . we will need to replace j3 with its Yosida rcgulariza- 
tion j3 T : M — > R, namely the nondecreasing, Lipschitz continuous derivative of the convex C 1 function 
(3 T (x) := min ye R{|?/ — x\ 2 /2t + /3(y)}, cf. e.g. [3J[TU]. In Problem 13.61 below, we set the regularization 
parameter equal to the time-step, in view of passing to the limit simultaneously in the time discretization 
and in the Yosida regularization as r — > 0. Furthermore, we will have to work with a suitable truncation 
of the coefficient a(X) in f|2 .58[) . cf. Remark 13 . 71 below for further comments. 

Problem 3.6 (Time discretization of the irreversible isothermal system). Starting from the triple of 
data (u^u^ 1 ,^) defined as in (|3.2[) and considering the discrete approximations (Q* T )?J_i °f the given 
temperature profile 6*, find {u^, X k , ( k }f=i G H 2 (n;R d ) x x L 2 (n), such that for all k = 

1 , . . . , K T there holds 

+ V ( K [{a(X k T ))+ + S p—^L + £ (b((X k )+)u k ) = f k a.c. in 0, (3.21) 

^—2^1 + g + s (x *) + ^(x*) + 7(x fc) = -^((x^- 1 )+) g(u "' 1) y (u "" 1) + e*^- 1 a.c. in n, 



(3.22) 

/ X k — "X k ~ 1 \ 

Cr ea( T — a.e. in n. (3.23) 



Remark 3.7. In Problcm l3.6l we need to approximate (3 by a Lipschitz continuous function j3 T because, 
only with such a regularization can we test equation (|3.22[) by the discrete difference T~ 1 (( r B(X k ) + 
(3 T (X k )) — (23(X^ _1 ) + f3 T (X k ~ v ))) (cf. the following Seventh a priori estimate). Hence, we need to take 
the positive part of a in (|3.21[) because, replacing f3 by its Lipschitz regularization f3 T , we are no longer 
able to enforce the constraint that X k £ [0, 1] a.e. in Q. Therefore, at the discrete level we loose all 
positivity information on the coefficient a(X). The lack of the constraint X k € [0,1] also motivates the 
truncations b((X k ) + ) in (|3.21|) and b'((X k ) + ) in (|3.22l) . mainly due to technical reasons (cf. the First a 
priori estimate). 

3.2 Existence for the time-discrete problems 

First, we prove the existence of solutions to the semi-implicit schemes (|3.17p - (|3.19[) and (|3.21[) - (|3.23D . 

Lemma 3.8 (Existence for the time-discrete Problems 13.41 and 13.61 ^1 = 1, p = 0). Let \i = 1 and p = 0. 
Assume Hypotheses (I)-(V), and (|2.47[) - (|2.51[) on the data f, g, "do, Uo, Vo, Xo- 

Then, Problems \3.4\ and [3~d\ admit at least one solution {(w k , u k , X k , Cr)}k=i ant ^ {( u t> ^1 (r)}k=i> 
resp. 

Furthermore, if (|2.84[) holds, then any solution {(w k , u k , X k , Cr)}£i of Problem \3.4\ fulfills 

w k (x) > w = hi^o) > for a.a.xe CI. (3.24) 

Proof. We treat Problems 13. 4[ and 13.61 in a unified way, and proceed by induction on k. Thus, starting 
from a quadruple (u k ~ 2 , w^ 1 , u*" 1 , X^ 1 ) G H$(fl;R d ) x ff^fi) x H 2 {fl;R d ) x W hp (n), we show that 
there exist functions (w k ,u k ,X k ,( k ) and (u k ,X k ,( k ), resp., solving (|3.17j) - (|3.19|) for Problem IIOI and 
(|3T2I jl -(|3T25 |l for Problem 1531 resp. 
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Step 1: discrete equation for X. In the irreversible isothermal case (i.e. for Problem 13. 6p . in order 
to solve (|3.22p we start from the approximate equation 



-yk -yk—1 / -yk -yk — 1 

A r,e a t / A i\e — a t 



T 



/ (3.25) 
= _ 6 , (x ,_ 1)£ (u^)^ £ (u^) + Q ^ a e . n ^ 

where e > and a e is the Yosida regularization of the operator a. Clearly, p.25[) is the Euler equation 
for the minimum problem 



mm 



where the function 



X - X 



fe-i 



2 



fc-1 



a £ ( X X _ T ) ] da- + $(X) + / (/3 r (X)+7(X))da;+ / /i^Xdz 

n 



h k - x := 6'(X^- 1 )e(u^- 1 )R e e(u^- 1 )/2 - 0**" 1 is in L 2 {9). (3.26) 

The latter admits a solution X^ e by the direct method of the calculus of variations (also taking into 
account the fact that f3 T is bounded from below because f3 is). We now want to pass to the limit in 
(I3.25|) as e I 0. Note that, a comparison in (|3.25p and the fact that a e is Lipschitz continuous yield 
that B(X* e ) e L 2 {Q). Then, following Sec. 3] (to which we refer for all details), we multiply (j3~23|) 
firstly by X\ e — X^ _1 , and secondly by CB(X^ e ) — 23(X^ _1 ). To perform the latter estimate, we rely on 
the Lipschitz continuity of /3 T and 7, as well as on the monotonicity of a e , yielding 



/fe-i 



a 



Xt - £ Xt ] (s(x* s ) - stxjr 1 )) dx > 0. 



It follows from these tests that there exists a constant C > 0, depending on r > but not on e > 0, such 
that 

sup(||x^|U 1 , P(0) + ||a3(x^)|| L2(0) )<c. 

£>0 

By comparison, sup e>0 ||a £ ((x£ e — X k .~ 1 )/r)\\ L 2^ < C. Also in view of the regularity result (|2.28[) . there 
exist (X*,Cr) e W 1+CT ' P (Q) x L 2 (fi) for all < ct < 1/p such that, up to a subsequence, (X* e ) e strongly 
converges in W 1,v {Vt) to X\ as e — ► 0, and (a £ ((Xji e — X^ _1 )/r)) e weakly converges in L 2 (S1) to £™ as 
£ — >• 0. Therefore, 



£^0 Jo 



so that Cr € "((^r - X r _1 )/ T ) thanks to [3J p. 42]. Thus, passing to the limit as e -> in (|3T2"5)) for r > 
fixed, we conclude that the functions (X*,£*) fulfill (|3.22|> . 

Clearly, the direct method of the calculus of variations also yields the existence of a solution to the 
minimum problem (13.191) . 



Step 2: discrete equation for u. Next, we solve (|3.18[) . which can be rewritten a.e. in f2 in the form 

(Id + tV ((a(X k ) + 5)-)+ r 2 £ (6(X*) •)) (uj) = r 2 f* + rV ((a(X*) + ^u^ 1 ) + 2a*" 1 - u^ 2 . (3.27) 

Combining the fact that X^ £ [0, 1] a.e. in with (|2.18[) on a and b and (|2.3[) - (|2.4p . we conclude that (the 
bilinear form associated with) the operator on the left-hand side of the above equation is continuous and 
coercive. Hence, by Lax-Milgram's theorem, equation Q3.27P admits a (unique) solution u£ £ HQ(fl;M. d ). 
Since the right-hand side of (|3.27p is in L 2 (il;R d ), relying on the regularity results of, e.g., [12], we 
conclude that in fact u£ £ i?Q(f2;R d ). The analysis of (|3.2ip follows the very same lines. 
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Step 3: discrete equation for w. Finally, let us consider the functional Q k 1 : 

s5r» : =J- / \w-w k ;- 1 \ 2 dx+ f e^- 1 ) f x "~_ x " - ) dx 



Q JO 
1 /" - ..k-1 



Now, Sr _1 is lower semicontinuous w.r.t. the topology of L 2 (Vl). Furthermore, in view of (|2.45[) and of 
the Young inequality we have for a fixed g > 

9*-» >-^|Mli 3(n) + yllV^||i 2( o) - elkll^ ( o) (3 2g) 

- C^U-Yma) + U\\% H ay + II (X* - X^/rU^). 
Choosing g sufficiently small, we thus obtain that there exist two positive constants c and C such that 

>c|Ml!ri(n) - C(l + ll^-'lll^n) + \M H ay + IK** - ^Mh^ for all w E H^Q) . 

This shows that the sublevels of 9t _1 are bounded in if 1 (17). Hence, again by the direct method in the 
calculus of variations, we conclude that there exists w k E ATgmin w&H i/^S^ 1 (w), and satisfies the 
associated Euler equation, namely p. 171) . 

Step 4: positivity. Let us assume in addition that (|2.84j) holds, and prove (|3.24[) by induction on k. 
Preliminarily, we prove by induction on k that 

w*{x) > for a.a. x E 17 and for all k E N. (3.29) 

Clearly (|3.29p holds for k = thanks to (|3.15p . It remains to show that, if w k ~ x > a.e. in 17, then 
Wr > a.e. in 17. Indeed, let us test ()3.17|) by — (w k )~ . Taking into account the definition (|2.37p of 0, 
we have that f n &(w k ^ 1 )(X k - X k ~ 1 )(-(w k )~) dx > 0. Combining this with the inequality 

-/(«£- ^- 1 )(-(^)-)do: > f / (|(^)-| 2 - {(w^ri 2 )^, (3.30) 
T Jn 2r J n 

and noting that (uj^ 1 )" = a.e. in 17, also in view of (|2.45[) we obtain 

-L f \(w L T )-\ 2 dx + c 2 f \W(w k T )-\ 2 dx<- f g k T (w k T )-dx<0, 
^ T Jn Jn Jn 

yielding (w k )~ = a.e. in 0, whence p.29|) . 

Now, to prove (|3.24[) . we observe that (|3.24p holds for k — due to (|3.15p . Suppose now that w k_1 > w 
a.e. in 17: in order to prove that w k > w a.e. in 17, we test p.!7p by — (w k — Wq)~- With analogous 
calculations as above we obtain 

^ / \(w k ~w )-\ 2 dx + c 2 j |V(^'- mo )-| 2 dx<- j g k {w k - m )-dx 
ZT Jn Jn Jn 

+ I e(^)(x^ 1 -x^(-(^-^ )-)dx<o, 

Jn 

where the last inequality is due to the fact that g k > a.e. in 17, and that 9(^)(X^" 1 - X*) > a.e. m 
17 by the previously proved (|3.29p and the irreversibility constraint. Thus, we conclude (|3.24p . □ 

The existence result for Problems 13.21 and 13.31 reads: 

Lemma 3.9 (Existence for the time-discrete Problems 13.21 l3~3l fi = 0). Let /x = 0. Assume Hypotheses 
(I) and (III)-(V), and (|2.47p ~ (|2.5ip on the data f, g, i?o 3 u o> v o> X . Furthermore, 
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1. if p = 0, assume Hypothesis (II); 

2. if p ^ 0, assume Hypothesis (VIII) and in addition that wq G L 2 (fl). 
Then, Problem \3.2\ admits at least one solution {{w k , u k , X k , ^r)}k=i • 

Moreover, ifg>Oa.e. infix (0,T), and Wq(x) > for a. a. x G fl, then any solution {(w k ,u k ,X k )}^ : l 1 
of Problem\KE fulfills 

w k {x)>0 fora.a.xeCl. (3.31) 

Proof. Step 1: existence of solutions. Our argument relies on existence results for elliptic systems 
from the theory of pseudo- monotone operators which can be found, e.g., in [3H1 Chap. II]. Indeed, we 
observe that system (|3.3|) — (|3.5f) can be recast as 



w k + (X k T - X^Mw*) + tA^-! (w l t ) = w k T - x + rg k in H^fl)', 

u k T + tV ((a(X k ) + 8)(u k - u^ 1 )) + r 2 £ (b{X k )u k ) = 2u k ^ - u^ 2 + r 2 f T fc a.e. in fl, (3 32) 

X k + rS(X^) + Tf3(X k ) + T1 {X k ) - T&(w k ) 3 X^ 1 - Tb'{X k - l ) e -^ ^fei^I % 1 a . c . in fi. 

Denoting by Rk-i the operator acting on the unknown (w k ,u k , X k ) and by Hk-i the vector of the terms 
on the r.h.s. of the above equations, we can reformulate system (|3.32[) in the abstract form 

Sn^.u^^iJ^. (3.33) 

In fact, mimicking for example the calculations in [471 Lemma 7.4], it can be checked that Olk-i is a 
pseudo-monotone operator (according to gSJ Chap. II, Dcf. 2.1]) on H 1 ^) x H^(Cl;R d ) x 
coercive on that space. Therefore, the Leray-Lions type existence result of [3H1 Chap. II, Thm. 2.6] 
applies, yielding the existence of a solution (w k ,u k ,X k ) to (I3.33|) . 

Step 2: non-negativity of w k . Let us assume in addition that g > a.e. in fi x (0,T) and wq > a.e. 
m 0. Then g k > a.e. in f2. To prove (|3.3ip . we proceed by induction on k and show that, if w k 1 > 
a.e. m n, then w k > a.e. in O. Indeed, let us test (|3.3[) by —(to*) . Taking into account the definition 
(|2"1T7)) of 8, we have that 



9«) ( (X* - X*"*) + pdiv ( "* "* 1 ] ) (-(w*)-)cLb = 



in \ \ ' / j 

(here we have kept p £ M also to encompass the case with thermal expansion, cf. below) . Combining this 
with the inequality (|3.30[) and noting that (w^ 1 )^ = a.c. in 0, also in view of (|2.45[) we obtain 



-L / |(^)-| 2 d.T + c 2 f |V(^)-| 2 dx<- [ g k (w k )-dx 
^ T Jn Jn Jn 



< 0, 



yielding (w k )~ = a.e. in ft, whence p.31[) . Under the additional Hypothesis (VIII) (which gives p.9[l ). 
an analogous proof of existence of solutions can be given for Problem 13.31 hence we omit to give the 
details. □ 



3.3 A priori estimates 



Notation and auxiliary results. Hereafter, for a given Banach space X and a AVtuple (b*)^ C X, 
we shall use the short-hand notation 

i k _ i k-1 uk _ Ofcfc-l 4. hk-2 

D Tlfc (6) := ^ T —, D 2 . fe (6) := D T)Jk (D T , fc (6)) = ~ ^ T , T ■ 

We recall the well-known discrete by-part integration formula 

J2 TT>r,kQ>>r = b^v 1 ^ - b° T vl - J2 rb^D^iv) for all {b k }^Z v K}f=i C X. (3.34) 

k=l k=2 
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Wc consider the left-continuous and right-continuous piecewise constant, and the piecewise linear in- 
terpolants of the values namely the functions 



b T :(0 1 T)^X defined by b T (t) := b*, 
b T :(0 1 T)^X defined by b r (t) := 6*-\ 
b T : (0,T) X defined by b T (t) := 



We also introduce the piecewise linear interpolant of the values {(&£ — b\ (namely, the values 

taken by the -piecewise constant- function b' T ), viz. 



b T :(0,T)^X b T (t):= 



t-t 



k—luk lM — 1 



tt-tb 



k-l _ uk-2 



forte {t k ;~\tp 



T 



Note that b' T (t) = fc (6) for t € (i* 

In view of (|2.47|) . (|2.48[) . and (|2.88|) . it is easy to check that the piecewise constant interpolants (fr^Zii 
(Sr)£=i> and (0*r)fcii of the values f r . 9r (EH), and 9*^ fulfill as r | 



f T ->f 

||^t@*r|U 1 (0,T;i2(n)) < 2||9 t 6*|| L i( ,T;L 2 (f2)) 



in L 2 (0,T;L 2 (n-,m. d )), (3.35a) 

in L 1 (0,T;£ 1 (fi))nZ 2 (0,:T;.H' 1 (fi)'), (3.35b) 

in 1/(0, T; L 2 (il)) for all 1 < p < oo, (3.35c) 

for all r > 0. (3.35d) 



Finally, we shall denote by t T and by t T the left-continuous and right-continuous piecewise constant 
interpolants associated with the partition, i.e. t T (t) := t\ if t\~ l < t < t k and t T (t) := t^T 1 if t^ 1 < t < 
t*. Clearly, for every t g [0, T] we have t r (t) 1 t and t r (t) 1 i as r -> 0. 

Propositions 13.101 and 13.111 collect in the cases p = and p ^ several a priori estimates on the 
approximate solutions, obtained by interpolation of the discrete solutions to Problems 13.21 I3T41 13.61 and 
Problem 13. 31 respectively. 

Proposition 3.10 (p e {0, 1}, p = 0). Let p = 0. Assume Hypotheses (I)-(V) and (|2.47[) - (|2.5ip on the 

data f, 5, $0j u o, vo, Xo- Then, 

1. in the case p G {0, 1} t/iere ezjst a constant S > smc/i £/ia£ /or t/ie interpolants of the solutions to 
Problem \3.2\ and to Problem \3.4\ there holds: 

(3.36) 
(3.37) 
(3.38) 
(3.39) 
(3.40) 
(3.41) 

(3.42) 
(3.43) 

(3.44) 



sup 

T>0 


|u T 


\ H 1 (0,T;H? 1 (£l;B. d ))nW 1 ^ (0,T;Ho(^ d )) — ^' 




sup 

T>0 


|u T 


I.L° o (0,T;flJ(fi;R< i )) < 5, 




sup 

T>0 


|u T 


lff 1 (0,T;L 2 (n;R d )) < 5, 




sup 

T>0 


\X T 


U a °(0,T;W I 'i'(n)) — S> 




sup 

r>0 


\X T 


L°°(0,T;W 1 '!'(O))nff 1 (0,T;L 2 (f2)) < «5, 




sup 

r>0 


\w T 


L°°(0,T;Li(^)) - ^' 




sup 

T>0 


\w T 


U'-(o,T;W 1 .'-(n)) < <5 for every 1 < r < 


- 2 

-1' 


sup 

r>0 


\w T 


BV([0,T];M/i.'''(a)*) - $ 




sup 

T>0 


|e(w T )|| L 2 +e ( 0iT;i 2 +e (n)) < S forany0<e< 


cr(d + 2) 

d 
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2. if p = in addition there exists S' > such that 



sup (||S(X T )|| i 2 (0 T;i 2 (n)) + ||CtIIl 2 (o.T:L 2 (o))) < S". (3.45) 

T>0 

Moreover, if <f> also fulfills Hypothesis (VI), then 

sup ||X T || L 2 (0 . T-w^+'-vin)) < <S" / or ever V < a <-. (3.46) 

r>0 P 

3. in the isothermal case with fx = 1, ifb" = (cf. (|2.86[) ) and (f> also fulfills Hypothesis (VI), estimates 
(|3.36[) - p.40p ZioZd. Moreover, there exists S" > swc/i that for (the interpolants of) the solutions 
to Problem^ 

sup (||S(^ r )|| i « (0iT . ia( n)) + ||/?r(X T )|| L oo (0)T . L2(fi)) ) < S", (3.47) 

T>0 

sup ||X T || L ~( 0i T;W 1 +-.p(n)) < S" for every < a < -, (3.48) 

T>0 P 

sup||X r || w ,i,<x J(0iT;L 2 (J2)) < 5", (3.49) 

T>0 

sup||C T |Uoo (0jT:i 2 (n)) <5". (3.50) 

T>0 

TTie constants in (|3.42[) . (|3.44[) . and (|3.46[) . (|3.48[) a/so depend on the parameters r, e, and a, respectively. 

Proposition 3.11 (p = 0, /) ^ 0). Let /i = and p ^ 0. Assume Hypotheses (I), (III)-(V), and 
Hypothesis (VIII); suppose that the data f, g, -do, Uo, Vo, Xo comply with (|2.47[) - (|2.5ip . and m addition 
that Wq € i 2 (J7). Then, for the interpolants of the solutions to Problem \3.3\ estimates (|3.39[) - (|3.4ip Zio/d 
wit/i a constant independent of M , whereas estimates (|3.36[) f|3.38[) . (|3 .45[) and (under the additional 
Hypothesis (VI)) (|3.46[) hold for a constant depending on M. Moreover, there exists a constant S = 
S"'(M) > such that 

sup ||w T || L 2 (0;T . ff i( O))nL oc (0:T;i 2 (o)) < S 1 , (3.51) 

T>0 

sup\\w T \\ H im T;H i(py) < S . (3.52) 

T>0 

We will treat the proofs of Propositions 13.101 and 13.111 in a unified way, developing a series of a priori 
estimates. 

Proof of Proposition 13.101 Most of the calculations below will be detailed on the discretization 
scheme (|3.3[> — (]3.5|) for the full reversible system, and whenever necessary we will outline the differences 
in comparison with the discrete systems of Problems 13.41 and 13.61 Furthermore, for each estimate we will 
specify the values of the parameters p and p for which it is valid and, to make the computations more 
readable, we will illustrate them first on the time-continuous level, i.e. referring to system (|2.57[) ~ (|2.59[) . 
First a priori estimate for p G {0, 1}, pel: we test ([238")) by u t ((2371) by 1, (12391 by X t , add them 
and integrate in time. This is the so-called energy estimate. We test (|3.4[) by u£ — ufjr 1 . Note that 

r f D* fc (u) ■ D T , fc (u)dx > i||D T)Jk (u)||| 3(n) - i||D Tlfc _i(u)||£, (n) (3.53) 

for all k = 1, ... , K T . Since X k T e [0, 1] a.e. in ft, thanks to (|2.18[) we have that a(X k ) > a.e. in fi, thus 
by (|2.3[) we have 

(V ((o(X*) + 6)D T , k (u)) , - u^ 1 )^ t (n) > CM^rM (u)||^ (n) . (3.54) 
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On the other hand, using that ||&(x£)||x,°o(n) < ||^||l°°(o.i) an d taking into account (|2.4[) . we find 
| (£ (fc(X*)u£) ,u* - u*" 1 )^ | < C 2 T||b|| i - (0) i)||uJ|| ff i ( n)||D T , fc (u)|| ff i (n) 



(3.55) 

^CiST\\D T , k (u)\\% Ha) + CT\\u° T f m{Q) + C S T\\u k ; - u° T f H1{Q) . 



We estimate the latter term by observing that 

A- 



„ = II^K-ur 1 )!!^^ <fcT 2 ^||D rj (u)||^ 1(0) <rr^||D TJ -(u)||^ 1(n) . (3.56) 



,k .,0||2 

1 r\\H 1 (n) 

3=1 3=1 3=1 



Altogether, collecting p.53[) - (|3.56[) and summing over the index k = 1, . . . , -frT r , we conclude 



-||D T ,K T (u)||i 2(0) + -C 1( 5^T||D^(u)||^ (J2) +p^T / 9(^)div(D T , fc (u))dx (3.57) 

j„ i i 1 ** fa 



fc=i fe=i 



< 2l|D r ,o(u)|| 

fe=l \j = l 

We multiply (|3.5[) by Xj: — X^T 1 . With standard convexity inequalities, we obtain 
r||D T , fe (X)||| 2(n) +$(%*)+ / ^)d.T+ / 7 (X*)(X*-XjT 1 )d a : 

< $(X^ J ) + jf ^(Xj-^da; +r jf D T , fe (X) (e(^) - ^'(X^Mu^R^u^ dx. 



(3.58) 



fc=i 

2 

fe=i 



We then test (|3.3p by r and add the resulting relation to (|3.57[) and (|3.58[) . summing over the index k = 
1, . . . , K T . The terms r J Q D T> k(X)Q(w k [) dx and pr J n Q(w k ) div(D Ti fc(u)) dx cancel out. Furthermore, 
we note that |6'(X^ _1 )| < C a.e. in ft since b G C 1 (M) and < X\ 1 < 1 a.e. in f2, and exploit the 
Lipschitz continuity of the function 7, which enables us to estimate the last term on the left-hand side of 
(|3.58[) . Ultimately, we obtain 

/ w^dx + ^r\\D T>k (X)\\l, {n) + ^{X^)+ f A?(X^)dz+i||D riifT (u)||| 2(Q) 

1 Kt 
fe=l 

< / w dx + $(X° T ) + f P{X T )dx + VtII^ii, 

+ ^Cr(||X^|j L2(n) + ||| £ (u^| 2 || L2(0) +l)||D T , fc (X)|| i2(0) 
fe=l 

A' T K T K T / k 

<C+-J2 T\\T>r,k(X)\\h(.a) +CY,r\\e(u k ;)\\t Hn) + c£> £ r||D TlJ -(X)||| a(n) 
fe=i fe=i fe=i \j=i 

+ c, EME T !i D ^( u )ii 2 ^( 

fe=l \3=1 

where the last inequality follows from assumptions (|2.48l) - (|2.5ip on the data, from the Young inequality, 
and from estimating t||X^||| 2 ,^ < 2r||X° H 2 ^^-, +2r||Xj: — X°|j 2 2 ^ and dealing with the latter term like 
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in (|3.56[) . Therefore, applying a discrete version of the Gronwall lemma (cf., e.g., [2H1 Prop. 2.2.1]), we 
conclude estimates (|3.39[) - (|3.41[) . as well as estimate 



sup ||u T || if i( 0jT . J ji( n . R d)) nM /i,oc(o,T;L 2 (a ; R< i )) < S, (3.59) 

T>0 

which in turn implies 

SUp||uT||i«(0,T;ffi(£);R<*)) < S. (3.60) 

T>0 

We can perform this energy estimate on Problem 13.41 as well: calculations (|3.53[) - (|3.57[) can be triv- 
ially adapted to (|3.18[) . whereas (|3.58|) derives from choosing in the minimum problem (|3.5[) X^T 1 as a 
competitor. We again conclude (|339]) - ([3"1I|) as well as ([339]) — ([3TBD]) . 

In the case of Problem l3.61 (|3.22p also features the term ( k , whence the additional term J n C^D T ^(X) dx 
on the left-hand side of (|3.58[) . Since £ a(0), by monotonicity the latter term in nonnegative. Taking 
into account this, replacing (3 with (3 T in Q3.58p . and observing that the coefficient of e(u k )R e e(u k ) on 
the right-hand side of p,19[) is bounded, we may repeat the same calculations as in the above lines. 
The coercivity estimate (|3.54[) goes through because a(X^), which is no longer guaranteed to be positive, 
is replaced by a(X k ) + . Furthermore, since X k < X^T 1 < X» < 1 a.e. in f2 (due to the irreversibility 
constraint), we have that (X k ) + £ [0, 1] a.e. in i7, thus we may again obtain (|3.55[) . 

Second a priori estimate for /j £ {0,1}, p = 0: following f8}j (see also \4^l ): we test (|2.58|) by 
— div(e(u t )) and integrate in time. We test (|3.4[) by — div(e(u£ — u^ -1 )). This gives rise to the following 
terms on the left-hand side: 

-t [ D? ifc (u).div(e(D Tifc (u)))dx>i / \s{T> T , k {u))\ 2 dx-\ / \e{T> T>k _ t {u))\ 2 dx, (3.61) 
Jn ^ Jn z Jn 

V((a(X*) + <S)D T , fe (u)) ■div(e(D T , fe (u)))dx = r f (S + a(X*)) div( £ (D T , fe (u))) ■ div(e(D Tlfc (u)))da: 



+ t f £(D T . fe (u))Va(X^) ■div(e(D Tifc (u)))da; 
Jn 

= /o+/i>«Cfr||D T)fc (u)||| 3(n) +/i, 

(3.62) 

the latter inequality due to (|2.9p . Moreover, always on the l.h.s. we have 

-r / £(6(X^K).div( £ (D r , fc (u)))dx = A 1 T /&(X*)A(D T , fe (u))-V(div(u*))d:z; 
Jn Jn 



+ 2A 2 r / 6(X^)div( £ (D r , fc (u)))-div(e(u^))da ; 
Jn 

+ Xit f dw(u k ;)\7b{X L T ) ■ A{D T , k {u))dx ( 3 - 63 ) 
Jn 

+ 2A 2 r f e(u k ;)Vb(X k T ) • div(e(D Tifc (u)))dx 
Jn 

= h + h + h + I 5 , 

(where A stands for the vectorial Laplace operator). On the right-hand side, we have 

-r ^f T ft -div( £ (D rife (u)))dz < C S T\\t k \\ 2 L2{n) + ^C|r||D T , fc (u)|| 2 ff2(n) , (3.64) 
where the latter inequality follows from (|2.9[) . We now move the integral terms I\, . . . , I5 to the right-hand 
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side. Let us fix < <r < 3 such that p > d + q (where p is the exponent in (|2.22[1 ). Then, 
|Ji|<r||div( e (D T , fe (u)))|| i = (n) || £ (D Tife (u))|| iti *- 5(n) ||Va(X^)|| ii+e(n) 

< jC|r||D T , fc (u)||^ (n) + 5r||e(D T , fe (u))||^_, (0) ||Va(X^)||i £i+ , (n) 

< ^ 3 2 r||D r , fc (u)||^ (0) + e 2 Cr||D T , fc (u)||^ (0) j|a'!| 2 L0O( _ m ^ m) ||VX^||i P(0) (3.65) 



+ C 2 d -C'||D rife (u)|| 2 c2(n) ||a'|| 2 Lo „ ( _ m . m) ||VX^||i P(0) 



< ^ 3 2 r||D T , fc (u)||^ (fi) 4 



where the first and second inequalities respectively follow from the Holder and Young inequalities, with 
d* as in (|2.12[) , the third one from (|2.13[) . and the last one taking into account estimates p.59[) for 
su Pk=i,...,K T |j D-7-,& ( u ) II i 2 (o) ; fl3-39[) for sup fc=1 m Kt ||X*||wi.i>(n)) which in particular yields that |x£| < m 
a.e. in Q x (0,T) for some m > 0, and from choosing g < C _1//2 (||a'||ioc(_ m m )S') _1 . Furthermore, taking 
into account that b(X%.) € L°°(0), one easily checks that 

\h+h\ < Cr||D Tife (u) ||ij2 (n) ||u^||H2 (n) 

A k (3 66) 

< -C 3 2 r||D T , fe (u)|| 2 H2(0) +C\\u° T \\ 2 H z in) +Cr5>||D TJ -(u)||^ (n) , 

where the second inequality follows from the Young inequality, from tHu^H^,,^ < 2t||u° + 
2r||uJ: — u° jj^aj-rj), and from estimating the latter term as in p.56[) . Analogously, again using that 
su Pfe=i k t ll-^TllvK 1 -p(n) < 5* and that |X*| < m a.e. in f2, we have 

\Ii + h\ < r ||D T ,*:(u)||/f2(Q) ||V6(Xj)|| £ s ( n)(|| div(u^|| i6(n) + \\e(u k ;)\\ L e {n) ) 

< C' r ll & 'llL~(-m,m)l|D Ti A ; (u)|| ff 2( n )||X T || l yl, P (0)||u T || H 2 (0 ) 



s 



(3.67) 



< gC3 2 r||D T , fc (u)||^ (n) + CS 2 ||u°|| 2 f2(0) + CS 2 t^t||D t ,,(u)|| 2 , 2(0) . 

3=1 

Collecting (|3.61|) - (|3.67[) and summing over the index k = 1, . . . , K T , we obtain 

i|| £ (D T ,^(u))||f, 3(n) + ^^r||D Tifc (u)||| 2(n) 

fc=i 

K T K T k 

< C + 2l| £ (D r , (u))|| 2 2(o) + C'Y, T\\f?\\i 2m + c£ r £ r||D Tj 

( u )IIh 2 (o)- 

fe=l k=l j = l 

Applying the discrete Gronwall Lemma once again, we conclude estimate Q3.36p , whence (|3.37[) . 

It is immediate to check that calculations (|3.6ip - (|3.67p can also be performed on the discrete momentum 
equation (|3.18[) in Problem 13.61 

Remark 3.12. The calculations for the Second a priori estimate carry over to the case the operator B 
is replaced by the s-Laplacian A s , provided that s > |. Indeed, this ensures the continuous embedding 
W s,2 (f2) C W 1 '^^) for some p > d, which is crucial in the above calculations, cf. (|3.65l) . 

Third a priori estimate for \i e {0, 1}, p = 0: BocCARDO&GALLOUET-iype estimate on (|2.57|l . As 
in the proof of [3H1 Prop. 4.2], we test equation (|3.3I) by H(w^.), where 

II : [0, +oo) -> [0, 1] is defined by U(w) = 1 - - — for some q > 0. (3.68) 

(1 + w) ? 
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Note that U(w k ) is well-defined, since w\ > a.e. in fi, and it belongs to n L°°(Q), as IT is 

Lipschitz continuous. Such a test function has been first proposed in [IB], as a simplification of the 
technique by Boccardo&Gallotjet [5]- We shall denote by IT the primitive of IT such that 11(0) = 
(hence II(u>) > for w > 0). Summing over k = 1, . . . , K T , we obtain 



< / n(^)d.T + ^r(|| 5r fc || il(f2) + ||D T , fe (x)eK fe )|| L1(n) )||n(^)|| L=0(n) , 
Ja k =i 

where the first inequality follows from (|2.45lh . the fact that VIT(w^) = ?(VwJ)/(l + w^Y +1 , and the 
second one from the convex analysis inequality J„ H(w^)(w k — _1 )do;> / n (II(to*) -n(w^ 1 ))da; and 
from the fact that, due to assumption (|2.49[) . we have 

n(w T )dx<c(\\w° T \\ Lim + i) <c. 

m 

Taking into account that < H(w^(x)) < 1 for almost all x E SI and all k = 0, . . . , K T , and relying on 
(j2~49j) and on (|3.35bl) . we conclude that 

E T / n EjC+i dT ~ g^r|lD T , fc (X)|| L2(n) ||e(^)|| L3( n) +C^. (3.69) 
fc=i Js:2 ^ + r ' 1 fc=i 

Now, wc argue in the very same way as in [JHl Proof of Prop. 4.2]. Combining the Holder and Gagliardo- 
Nirenbcrg inequalities (cf. (|2.11[) ]1 with the previously proved estimate (|3.41[) and with (|3.69[) . we see that 
(cf. Formula (4.35)]) 

VI < r < ^3C r , C; > OVr > : 5>||V«£||£ r(n) < C r £ T||D T , fc (X)|| ia(n) ||e(u£)|| i3(n) + C' r , 

k=l fe=l 

(3.70) 

where the restriction on the index r in fact derives from the application of the Gagliardo-Nirenberg 
inequality (|2.11j) . Next, for a sufficiently small e > such that a from (|2.16|l fulfills a > (2 + e)d/(d + 2), 
there holds 

||eK fe )||^ ( n) < c(ll«*llK& (n) + 1) < qi^ll^) (1 - 9)/ l«*llgt? ( Sr + c 

< CS^+^-^iS + \\Vw k ;\\ Lrm ) {2+eW ° + C (3.71) 
^ Mr, fenr „ ., d(d + 2) d+2 

(where wc have omitted to indicate the dependence of the constants on e and cr). The first inequality 
follows from (|2.44|) . the second one from the Gagliardo-Nirenberg inequality (|2.11j) with s = 2/ a and 
q = 1: in fact the constraints 

^>^. ^<'<£l 3 .e M =^-.(i-i) +1 -., <-) 

in accord with formula (|2.11l) . Finally, the last inequality in (|3.71[) is due to the Young inequality, 
with C e depending on the constant g > to be suitably specified, under the additional condition that 
r < (d + 2)/(d+l) fulfills 

- — < r. 



34 



Combining (|3.71[) with (|3.70[) . we immediately obtain 

E T ii Vw "n™ < ^E T !i D ^wiii 2 w + ^E T ii Vw "ii^) + c "- ( 3 - 73 ) 



K T K T K T 

T 

k=l k=l k=l 



Hence, we choose g > in such a way as to absorb the second term on the right-hand side into the 
left-hand side. Therefore, on account of p.40p sup T X)fe=i t II^ u; tII£' , (0) — which yields (|3.42p via 
(I3.4ip and the Poincare inequality. Finally, estimate p.44p ensues from p.42p and (|3.71[) . 

Observe that, when performing this estimate on the semi-implicit equation p,17[) . we will obtain on 
the r.h.s. of (EZD the term £f= x r\\ V^" 1 ||£ r(n) < 7-||Vt«or||^ (n) + Ef=x r||Vu;*||£ r(n) , and we can 
estimate r|| Vu^i-H^n) thanks to (|3.14j) . 

Fourth a priori estimate for p 6 {0, 1}, p = 0: comparison in (|2.58|) . It follows from estimates (|3.3G|) . 
p~3"Tj) . (12301), and from the regularity result (|2.6bl) . that 

SUP ||V ((a(X T ) + S)d t U T ) ||i2( ,T;£2( n;]R d)), Sup |j£ (b(X T )u T ) ||i~( 0l T;i=(n ; K<2)) < C. 
r r 

Thus, for p = estimate (|3.38[) follows from a comparison in (j3.4[) . 

The same argument carries over to p,18p and to (|3.21[) . 
Fifth a priori estimate for p £ {0,1}, p = 0: comparison in ()2.57p . In view of estimates and of 
p.35bp . a comparison argument in (j3.3[) yields estimate p.43p . The same for (|3.17[) . 
Sixth a priori estimate for p = 0, pgl: we test (|2.59p 6?/ 23(X) + /3(X) and integrate in time. We 
test p.5p by r (23 (X J;) + Arguing as for p.58p via convexity inequalities and referring to notation 
p.26p for the symbol ft* -1 , we get 

S(X*) + ^(x*)dB + r||S(X*) + £||| 8(n) 

<<&(X^ 1 )+ / ^- 1 )da ; + r|| 7 (X^) + ^- 1 || L2(n) ||a3(X^)+e T fc ||L 2 (0) 



< ^(X^ 1 ) + jf ftx^cte + ^r||3(X*) + £||| a(n) + Crdl^Hi,^) + 1) 

where the last inequality follows from < X^ < 1 a.e. in CI, and the fact that 7 is Lipschitz continuous 
on [0, 1]. Summing up the above inequality for k = 1, . . . , K T and taking into account a priori estimates 
P-37P and p.39p . we conclude 

sup||I5(X T ) +? T |U2 (0 , T;i 2 (n)) < C. (3.74) 

r>0 

From this bound, exploiting the monotonicity of ft and applying [101 Prop. 2.17], we deduce p. 4511 . In 
view of (|2.28p . from the estimate for 23 (X T ) we deduce p. 461) . 

Seventh a priori estimate for p = 1, b" = 0, and in the isothermal case: we test (I2.59P by 
d t {'B(X) + ft(X)). Since 6" = 0, we have that 6'(X^ 1 ) = b on Q. We test p~2"2)) by rD T)fc (2(X) +A-(X)) = 
(S(X*) + /9 r (X*) - 0B(x£ _1 ) + /^(X^ 1 ). We observe that 

/ 6 := / (X^-X^ 1 )(S(X^)-S(X^ 1 ))da;= f (VX^-VX*" 1 )-^, VX*)-d(z, VX^ 1 ))dx > 0, (3.75) 

and, if (|2.26[) holds, we have in addition 

h > c 7 / |V(X* - X*" 1 )^ = c 7 r [ r p - 1 |VD r , fc (X)| p dx (3.76) 
Jo Jn 

Moreover, by monotonicity we have 

/ (X k T - X k T - r ){ft T {X k T ) - ftAX^dx > 0, r f £(ftAX k T ) ~ ftAX^dx > 0. (3.77) 
Jn Jn 
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Furthermore, always by monotonicity, we get 

C*CB(X*)-3(X* _1 ))da; > 0. 



T 



it 



Here, in order to perform a rigorous argument we should approximate the graph a with a Lipschitz 
continuous function a e as we have done in Lemma 13.81 However we prefer not to do it now in order not 
to overburden the calculations. Combining (|3.75[) and (|3.77p with the inequalities 

(S(X*) + T (X*))(rD Tifc (S(X) + T (X)))dx 



J Sfi J Si 

and summing over the index k = 1,...,K T , we get (cf. (|3.26[) for the notation h^ 1 , with 9(u)^ _1 J 
replaced by ©*^ _1 ) 

\ f \n^)+P T {X^)\ 2 dx<\ [ \V(X°)+f3 T (X T )\ 2 dx 



]Tr / (ft*- 1 - 7 (xJ))D T . fc (B(X)+ / 8 T (X))dx. (3.78) 
fe=l Jn 



= /7 

Clearly, the first term on the right-hand side of (|3.78[) is bounded thanks to (|2.87[) . Applying the discrete 
integration by part formula (|3.34j) . we find 

h = (3(X^) +P T (X^))(h^-i + 7 (X^)) - (3(X°) +/3 T (X° r ))(^ + 7 (Xl)) 

1£l (3.79) 
- £ r(23(X^ 1 ) + /? r (X^ 1 ))(D T , fc _ 1 (/ l ) + D T)fc ( 7 (X))) . 

Now, by the Lipschitz continuity of 7 on [0, 1], we have ||DT,fe(7(^))||L 2 (f2) < C||D Ti /c(X)||x,2(n). Further- 
more, we find 

\\DT,k-i(h)\\*M < IP^fc-iCO-)!!^^ +Cr|&|||| £ (u^ 1 )| 2 - | e (u^ 2 )| 2 || i2(0) 

where the second inequality also follows from the fact that 6' is constant. Collecting (|3.78[) - (|3.79[) and 
the above inequalities, we thus infer 

i r Kt 

- / |S(xf-) +/3^(X^)| 2 d a : < C + ^^(||D T , fc+1 (X)|| i2(n) +^)||CB(X^) 4-/3^(^)11^^), 

where we set ~D t> k t +i(X) = 0. Then, estimate Q3.47P ensues via the discrete Gronwall Lemma, taking 
into account that 

k t 

5>(||D T , fe (X)|| i2(0) +^) <C (3.80) 
fe=i 

in view of (|3.35dl) . (|5T55]) . and (|3"H)|) . Ultimately, (pQg)) follows from ([3~4T)l and the regularity re- 
sult (|2~28)) . 

Eighth a priori estimate for /j, = 1 and in the isothermal case: comparison in (|2.59p . From a 
comparison argument in (|3.22p , we conclude that D Ti fc(X) + Cr is estimated in L°°(0, T; L 2 (f2)). Then, 
jSHSD and ([3~50| follow from the fact that J Q D Tjfe (X)C£ dx > 0. □ 
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Proof of Proposition [37TT1 The a priori bounds (|3.39|) - (|3.41l) follow from the calculations developed 
for the First a priori estimate, which also yields (|3.59[) and (|3.60|) for a constant independent of M > 0. 
The BocCARDO&GALLOUET-type Third estimate is replaced by the following 

Ninth a priori estimate for p = 0, p 7^ 0: test (|2.57p by w. We test (|3.1ip by ru>£. Summing over 
k = 1, . . . , K T and recalling (|3.7p ~ p.l0p we obtain 



1 

2 

< 
< 



yf T | 2 da; + ci ^r 
fe=i 



,,0 |2 



dx 



E 

k=l 
K 



Ww^dx 



g*w*dx 



(|D T)fc (X)| + H| div(D T , fc (u))|) \Q M (w*)\\w*\dx 



I ^0 1 



'dx 



v 



fe=i 



„fc||2 



li? 1 (a) 



aE- 

fe=i 



ll^llff^n)' + lpT,*(X)||£ a( n) + P 2 || div(D Tlfc (u))||£ a(n) 



(3.81) 

for a suitably small constant v > 0, where we have used that the terms Qm{Wt) are uniformly bounded 
(by a constant depending on M > 0). Hence, estimate (|3.5ip follows from using that wo £ L 2 (£l), taking 
into account the previously proved bounds p.40p and (|3.59p , and applying the discrete Gronwall Lemma. 
Estimate (|3.52p then ensues from a comparison in (|3.1ip . in view of the previously proved estimates. 

Finally, relying on (|3.5ip . we are able to perform the analogue of the Second a priori estimate on the 
momentum equation in the case p 7^ as well, as the following calculations show. 

Tenth a priori estimate for p = 0, p 7^ 0: test (|2 . 58|) by — div(e(u t )) and integrate in time. We test 
(|3.12p by — div(e(uji — u^ -1 )). Every term can be dealt with like in the Second estimate, in addition we 
need to estimate the term 



rp I V(e M (<))div(e(D r , fc (u)))da; 
n 



< apr||v(e M (^))||| 2(n) 



Pr,fc(u) 



2 



Choosing v sufficiently small in such a way as to absorb ||D Ti /-(u) 



into (|3.62p . and estimating 



|| \7(<d M (w^))\\ 2 L 2(m vm <|3-5ip (observe that 6m is Lipschitz continuous), we re-obtain (|3.36p - (|3.37D . for 



\L 2 ({1) 

a constant depending on M. Moreover, estimate f|3 .38[) ensues from a comparison in (|3 . 1 2[) . 

Finally, estimates p.45p and (|3.46p can be obtained by repeating on equation (|3 . 1 3[) the very same cal- 
culations developed for the Sixth estimate: again, we get bounds depending on the truncation parameter 
M. □ 



Remark 3.13. A close perusal of the proof of Proposition 13.101 and in particular of the calculations 
performed in the Second and Fourth a priori estimates, reveals that in fact estimates (|3.39p ~ p.43p hold 
for constants independent of 5 > in both cases p, = and p = 1. This will play a key role in the proof 
of Theorem [5] 

We conclude this section by mentioning in advance, for the reader's convenience, that the relevant 
estimates 

1. for the proof of Thm.[T]are the First, Second, Third, Fourth, Fifth, and Sixth a priori estimate; 

2. for the proof of Thm.[2]are the First, Fourth, Sixth, Ninth, and Tenth a priori estimate; 

3. for the proof of Thm.[¥]are the First, Second, Third, Fourth, and Fifth a priori estimate; 

4. for the proof of Thm. [S] are the First, Second, Fourth, Seventh, and the Eighth a priori estimate. 
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4 Proofs of Theorems [T], [2], and [3] 



4.1 Proof of Theorem CO 

Preliminarily, we rewrite equations p.3[) - p.5[) in terms of the interpolants w T , w_ T , u r , u r , u T , u T , X 7 
X T , X T , and £ T , namely 



/ / uJ r (p t dxds + / / dtX T Q(w T )ipdxds + / / i^(«; T )VMT T V</?da;ds 

Jo Jo Jo Jsi Jo Jn 

= [ [ g T tpdxds- [ w T {t)tp{t)dx + [ w Q ip{0)dx for all ip g t€[0,T], 

Jo Jn Jn Jn 



(4.1) 



d t u T (t) + V ((a(X T (t)) + S)d t u T (t)) + £ (b(X T (t))u T {t)) = f T (t) a.e. in f2, for a.a. t g (0, T), (4.2) 

a t x T (t) + s(x r (t)) + ? T (t) + 7 (x T (t)) = -b\x T (t)) e ^ {m f ( - Tm + e(w T (t)) 

z (4.3j 
a.e. in ft, for a.a. t g (0, T), 

where for later use in (|4.1| we have already integrated by parts in time, and & is as in (|2.57p . In what 
follows we will take the limit of (|4.1j) - (|4.3[) as r 4- by means of compactness arguments, combined with 
techniques from maximal monotone operator theory 

Step 1: compactness. First of all, we observe that due to estimates (|3.36[) and (|3.38[) . there holds 
||u T - u r || L 3o( 0iT;if 2( n . R d)) < T 1/ ' 2 ||9tU T || i 2( 0;T . Jf 2 (n . RC i )) < St 1 / 2 , 

||u r - d t U T \\ L ^^ !T ;L 2 (Q;R d )) < rl/,2 ||^Ur||L 2 (0,T ; i 2 (f!;R li )) < St 1/2 . 

Therefore, (|3.36[) — (|3.38|> . joint with (|4.4I) and well-known weak and strong compactness results (cf. [51]), 
yield that there exist a vanishing sequence of time-steps (r^) and u as in (|2.53[) such that as k — > oo 

u Tfc ^*u in i7 1 (0,T;7J^(rj;]R d )) n W 1,oo (0, T; H^(Cl; R d )), 



u rfc , u Tfc , u Tfc -> u in L°°(0,T;iJ 2 - e (fi;]R d )) for all e g (0,1], 
5 t u Tfc d tt n in L 2 (0, T; L 2 (£l; R d )), 



as well as 



(4.5) 

<9 t u Tfc -> <9 t u in L 2 (0,T;ff 1 (fi;l ,! )). 

A stability estimate analogous to the first of (|4.4[) . the a priori bounds p.39[) . (|3.40[) and (|3.45|1 and 
the previously mentioned compactness arguments, imply that there exist X g L°°(0,T;W 1 ' p (SY)) n 
i? 1 (0, T; i 2 (fi)) and (,A £ L 2 (0, T; L 2 (£7)) such that, along a not relabeled subsequence there hold 
as k —> oo 

X Tfc -*X in ^(0,^^(0)) R 1^(0, T;L 2 (fi)), 

X Tfc , X rfc , X Tfc in ^(O.Tjr-'*^)) for all ee (0,1], 

£ Tfr -£ inL 2 (0,T;L 2 (r!)), (4.7) 
!B(X T J^A in L 2 (0,T;L 2 (ft)). (4.8) 

Furthermore, if in addition <fi complies with (|2.26|) and (|2.27|) . then, due to (|3.46[) we also have the 
enhanced regularity (|2.63[) . and the strong convergence 

X Tk , X Tk , X Tk -> X in L s (0,T; for all 1 < s < oo. (4.9) 

As for (w T ) T , estimates (|3.41|l - (|3.43|l and a generalization of the Aubin-Lions theorem to the case of 
time derivatives as measures (see e.g. [48j Chap. 7, Cor. 7.9]) yield that there exists w as in (|2.52|) such 
that, up to the extraction of a further subsequence, as k — > oo there hold 

w Tk ,w Tk ^w ini7(0,T;W^(fl)), 

w Tk ,w Tk ^w in L r (0,T;W 1_e ' r (n)) n Z/(0, T; L X (Q)) for all e g (0,1] and 1 < s < oo. l ' J 
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Furthermore, by an infinite-dimensional version of Holly's selection principle (cf. e.g. g]) we have w Tk (t) — i 
w{t) in W iy (n)* for all t £ [0,T]. Taking into account the a priori bound ([3~4Tj) of {w Tk (t)) Tk in L^fi), 
we then conclude that 

w Tfc (t) in M(O) for all i G [0,T]. (4.11) 

Clearly, the second of (|4.10l) implies that w Tk —> w a.c. in il x (0, T), hence by the continuity of O we 
also have 9(w Tfc ) — > 9(«;) a.e. in x (0,T). Moreover, estimate (|3.44[) guarantees that (Q(w Tk )) Tk is 
uniformly integrable in L 2 (0, T; L 2 (f2)). Therefore, thanks e.g. to [121 Thm. III. 3. 6], we conclude that 

e(w Th )^Q(w) mL 2 (0,T;L 2 (n)). (4.12) 

Step 2: passage to the limit in ([4TT] ) -(|4.3 [) . It follows from (|4.6|l and (|4.12|) that 

d t X Tk Q(W Tk ) ^ X t e(w) inL^T^fi)). (4.13) 

Moreover, P~TUj) and ([2~T71) easily yield that (K(w Tk )Vw Tk ) Tk is bounded in L r (0, T; Z r (0)). Since 

K{w Tk )^K(w) in i s (0,T;L s (fi)) for every s G [l,oo) (4.14) 

taking into account (]4. lOf) . we can pass to the limit as k — > oo in the third integral in the first line of 
(|4.1[) . Convergences (07TDJ) (gHH), (|4TT3"|) - (|P4"|) , as well as (|3.35bl) for {g Th ) Tk , allow us to take the limit 
of (|4.1[) as Tfc | 0. Hence we conclude that (w, X) comply with (|2.57|) . 

As for the passage to the limit in (|4.2D . we observe that f)4.6[) . the compact embedding (|2.14p . and 
(|2~T8l) imply that a(X T J -* a(X) and 5(X Tfc ) 6(X) in L°°(0, T; L°°(f2)). Therefore, by the (03]) we im- 
mediately conclude that V ((a(X Tk (i)) + <5)<9 t u Tfc (t)) V {{a{X(t)) + S)d f u(t)) and £ (6(X Tfe (t))u Tfc (t)) 
£ (6(X(i))u(i)) in L 2 (tt) for almost all i G (0,T). Also relying on the third of flU) and on p.35aj) for 
(f Tfe )fe, we find that (u,X) fulfill ([235)l . 

Finally, combining (|4.6[) with (|4.7[) ~ (|4.8[) and taking into account the strong-weak closedness of the 
graphs of the operators <9/3, 23 : L 2 (il) =4 L 2 (f2), we immediately conclude that f|4.7[) — (|4.8[) hold with 

£(x, t) e P(X(x, t)) and X(x, t) = i)) for a.a. (x, t) G x (0, T). 

We also observe that (|43|) , (g^J), and (|2~Tg|) yield 

5 , ( (f))£ (u Tfc (t))R e£ (u Tfc (t)) ^ g (u(t))R eg (u(t)) . n L 2 (QiT . L 2 m (415) 

fc 2 2 

Therefore, relying on (|4.6[) and the Lipschitz continuity (|2.20[) of 7 on bounded intervals, we easily pass 
to the limit in (O and infer that (X,£) fulfill (|2TM| - ([2TCT)) . 

Step 3: proof of total energy equality (|2.62[) . We choose 93 = 1 in (|2.57[) . test (|2.58[) by u t and integrate on 
time, and test (|2.61[) by X t and integrate on time, then add the resulting relations. Some terms cancel out, 
and to conclude (|2.62l) we use the chain rule formula (|2.83l) for the functional ${X) :~ $(X) + J„ Vy(X) dec, 
as well as the fact that 

Qa ol (&(X)u,u)^ = ia c i(6'(X)X t u,u)+a ol (6(X)u,u t ) 

= lj b , m e(u)B^e(u) ^ + (£ for a a t £ % T) _ 

□ 



3!) 



4.2 Proof of Theorem |2j 

Outline. First of all, relying on the a priori estimates of Prop. 13.111 we pass to the limit as t — > in the 
time-discretization scheme (|3.11[) - (|3.13[) : we thus obtain the existence of a triple (%,Um,Xjk) solving 
the truncated version of system (|2.57p - (|2.59[) . 

Secondly, we perform the passage to the limit as the truncation parameter M tends to +oo. In this 
step, we need to obtain for the functions (wm)m a bound in the spaces specified in (|2.67[) . independent of 
the parameter M. In this direction, the key estimate consists in testing (|4.17[) below by (a truncation of) 
W]\j £ which is now an admissible test function for f|4. 17[) : it is indeed in view of performing this 

test, that we need to keep the two passages to the limit as r — > and as M — > oo distinct. In order to 
carry out the calculations related to such an estimate, we need to carefully tailor to the present truncated 
setting the formal computations outlined in Remark 12.101 

Step 1: passage to the limit ss t -> 0, for M > fixed, in (|3.11[) - (|3.13p . The argument follows the very 
same lines as the one in the proof of Thm. [TJ it is even easier, due to the truncations of the functions K 
and 0. Therefore, we omit the details. Let us just observe that passing to the limit as r — s- in (|3.11l) 
leads to a solution 

w M € L 2 (0,T;H 1 (fi))nC ([0,r];L 2 (fi))ni7 1 (0,T;i? 1 (fi)') (4.16) 
of the truncated enthalpy equation 

(w t ,(p) H i ( m+ / X t Q M (w)(pdx + p / div(u t )Q M (w)ipdx + / K M (w)\7w\7ipdx = (g,f) H i,m (4.17) 
Jn Jn Jn 

for all tp £ i? 1 (il). Indeed, regularity (|4.16[) follows from estimates (|3.51[) - (|3.52j) . and in turn it allows 
for the stronger formulation (|4.17p (in comparison with (|2.57|l ) of the enthalpy equation. Therefore, for 
every M > the (Cauchy problem for the) truncated version of system (|2.57[) - (|2.59p . consisting of (|4.17l) 
and of (|2.58p - (|2.59p with 9 replaced by ®m, admits a solution (wm, Ujf) ^m)i with the regularity (|4.16p 
for w M and (j2~551) - (j2~H|) for {u M ,X M ), further fulfilling ([2~B0"j) - ([2'3Tj) (with 9 M in place of 6). 
Step 2: passage to the limit as M — > oo. Let (wm, vlm, Xm)m be the family of solutions constructed in 
the previous step. Since estimates p.39p - (|3.4ip and (|3.59p - (|3.60p hold with a constant independent of 
M, we conclude by lower semicontinuity that 



3C>0 MM > : \\wM\\ L ^^ 0tT . L i^ + \\u M \\H 1 (o,T-Ml l (n-m d ))nw 1 ^(o,T-L 2 (a-,M d )) 

+ \\ x M\\L°°(o,T;W 1 -p(n))nH 1 (o,T-L 2 (n)) < C- 

Next, introduce the truncation operator 

!-M if r < -M, 
r if |r| < M, 
M if r > M, 

and the sets 

f «4f := {(x,t) eSlx(0,T) : \w M (x,t)\ < A/}, stf M := {x £ Q : (x,t) £ #/ M ] 
1 G M ■= {(x,t) £ n x (0,T) : \w M (x,t)\ > M} 7 G M := {x £ il : (x,t) £ ff M } ■ 



(4.18) 



(4.19) 



Hence, we test (|4.17p by 7m(wm) and integrate on (0,i), t £ (0,T): observing that 

K m (w m )Vw m -V(7m(-wm)) = K(7 m (w m ))\V('T m (wm))\ 2 1 x f . f . 

e M (%) = e(T M (^Af)) J ' ' ' h [ ' ) 
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wc thus obtain 



- / \7 M (w M (t))\ dx+ / / K{7 M (w M ))M7 M (w M ))\ z dxda 

z Jn Jo Jn ^ 

< I I \7 M (w M (0))\ 2 dx + f f \g\\T M [w M )\dx + f f \e M \\Q(?M(w u ))\\7 M (w M )\teda, 

Jo Jn Jo Jn 



2 



o Jn 



where we have used the place-holder £m '■= dt^M + pdiv(<9 t ujv/). Now, arguing in the very same way as 
in Rmk. 12.101 we observe that 

K(7 M (w M ) ) I V (T M (wm ))\ 2 dx >cj ( 1 1 VT M (w M )\\h m + \ I ?m (w M ) 1 1 S^, (n) ) ds - C 

(4.22) 

for positive constants c and C independent of M. Let us now focus on the on the r.h.s. of (|4.21[) : note 
that ||9"m(wm(0))|| 2 2(q) ^ II u; o|||2(q-), whereas the second integral term can be estimated thanks to (|2.48[) 
on g. Taking into account the growth (|2.43[) of and (|2.46[) . the third integral can be estimated by 

cf f \hi\(\TM(wM)\ q+1 +l)dxds< g [ f \7 M (w M )\ 2 ^ds + C g , (4.23) 
Jo Jn Jo Jn 

where we have used estimates (|4. 18[) . Choosing g > sufficiently small, we can absorb the integral term 
on the r.h.s. of (|4.23p into the r.h.s. of (|4.22j) . As in Rmk. 12.101 we thus conclude that 

3C>0VA/>0 : ||T M (wA/)|| i 2( T;ff i (n))nL oc( 0T . L 2 (n))nL 2(< ! +i) (0jT . L 6(< ! +i)( )) < (7. (4.24) 

We now use (|4.24[) in order to infer an analogous estimate for the family (u>m)m- To do so, we 
preliminarily observe that from the bound for ||Tm(wm)||l°°(o.t ; l 2 (q)) we infer 

C> [ \7 M (w M )(t)\ 2 dx > f M 2 dx = M 2 \0 t M \. (4.25) 
Jn J ei, 

Therefore, upon testing (|4.17[) by wm, integrating in time, and repeating the same calculations as above 
(also relying on (|3.9[0 . we end up with 



^ / \w M {t)\ 2 dx + c 10 [ ||Vw M ||£2 (n) ds + c ^ WwMf^t+l^^ds 

r-l 



(4.26) 



<C+- \w Q \ 2 dx + Q X I \\wM\\ 2 H i,mds + C ei I \\g\\ 2 m ds + h 
z Jn Jo Jo 

for some g\ > to be specified later, and we estimate 

\£M\\@M{wM)\\w M \dxds + / / \lM\\®M(w M )\\wM\dxds 

Jo Je* M 

M&V^d* ( 4 - 27 ) 



+ C e* J q \\ i M\\mn)\\®M(w M )\\ L 3 { ffs M) ds + g 3 \\w M \\ H i {&li) ds 

where we have argued along the same lines as in Rmk. 12.101 Now, observe that for almost all t £ (0, T) 

jAf| 2 / CT + l 



\®m(wm)\\ 2 lW < C(\Q(M)\ + |9(-M)|) 2 |^| 2 / 3 < C(\M\ 2 ^ + l)^ 2 / 3 < C 



Af 4 /3 



where the second inequality is due to the growth (|2.43p of 6, and the last one to (|4.25[) . Observe that 
2/cr - 4/3 < for d = 3 thanks to ([2~ToT) . therefore (\M\ 2 / a + 1)/M^ 3 ^ as M oo. For d = 2, 
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(4.30) 



in (|4.27|) taking into account the Sobolev embedding H 1 ^^) C L s {ff M ) for every s G [l,oo), we can 
replace ||0Af (wa/ )|| 2 3(^>s ) with ||0m(u'm)||^2+5(^s ) for any e > 0, and tunc e in such a way that the 
latter term will again converge to as M — > oo. Therefore the third term on the r.h.s. of (|4.27|) is 
bounded. It remains to choose g 2 in such a way as to absorb the term J Q * IIwmII^J+iv^s , d.s into the 
l.h.s. of (|4.26p . and gi, q% so that Q\ + is sufficiently small. Also applying the Gronwall Lemma, we 
conclude that 

3C*>0VAf>0: \\w AI \\ L 2 {0tT . H i m)nL ^ {0 ^ T . L 2 m) <C . (4.28) 
With easy calculations we also find 

3 C > VM > : \\w M \\ L 2 {q +i) {0 , T . L 6( q +i) + \\dtw M \\mo,T:W 2 -s(ny) < C, (4.29) 

with s as in (|2.67p . the estimate for dtWM following from a comparison in (|4.17[) . 
We are now in the position to obtain the further estimates 

3C>0 VM > : ||uA/||i/i(o,T;ff2 ( n ;R <i)) nW /i,oo(o,T;ff 1 (0;R d ))ni/ 2 (0,T;L 2 (f2;] 

+ ||23(XAf)||L 2 (0,T;L 2 (J2)) + 1 1 £m 1 1 1,2(0^1,2(0)) < C , 

where £ M is the selection in @(X M ) fulfilling ([2~6Tj) . Indeed, (|4~30)) can be proved relying on the previously 
obtained (|4.18|) and (|4.28j) - ([4.29p . by performing on the truncated version of system (|2.57|l - (|2.59|) . the 
time-continuous analogues of the Sixth and Tenth a priori estimates (cf. Sec. I3.3|) . 

Hence, we can carry out the passage to the limit argument as M — > oo. Relying on the bounds (|4.18[) 
and (|4.30[) and on the compactness tools already exploited in the proof of Thm. [TJ we find that there exist 
(w, u, X) and a (not relabeled) subsequence of (wm, umj Xm)m such that (the time-continuous analogues 
of) convergences (|4.5p ~ (|4.9p hold as M — > oo. Furthermore, estimates (|4.28p - (|4.29p . the aforementioned 
generalization [351 Chap. 7, Cor. 7.9] of the Aubin-Lions theorem to the case of BV-functions applies, 
and convergences (|4.10D improve to 

w M ^*w in ^(O.TjfT 1 ^)) n L^" +1 ^0,T; L a ^ +l \Q)) nL°°(0,T;L 2 (Q)), 

w M -> w in L 2 (0, T; W 1 ^' 2 ^)) n i CT (0, T; L 2 {il)) for all e € (0, 1] and 1 < a < oo. 

with w G BV([0, T]; VK 2 ' 5 (rZ)') and s as in ([2J7)l . Relying on (|4~3Tj) and on the Lipschitz continuity of 6, it 
is not difficult to infer that V8m(wm) — 1 V9(u)) in L 2 (0, T; L 2 (f2)) as fc — >• oo. Furthermore, combining 
(gUD and the last of (g3J| we also have that dW(d t u M )®M(w M ) div(u t )0(iu) in L 1 (0,T; L 3 / 2 (fl)). 
Therefore, we are able to pass to the limit in (|4.17JI (with test functions as in the statement of Thm. [2]), 
and in the corresponding equations for u and X in the case p ^ 0, which concludes the proof. □ 



5 Proof of Theorem [3] 

Let (ui,Xi), i = 1,2, be two solution pairs like in the statement of Theorem [3] and set (u, X) := (ui — 
u 2 ,Xi — X 2 ,). Taking into account that a is constant (cf. (|2.72p ). hence a(Xj) = a > for i = 1, 2 (cf. 
(|2~72l) s ). it is immediate to check that (u, X) fulfill a.e. in fi x (0, T) 

u« + £(&(Xi)u) + £((6(Xi) - 6(X 2 ))u 2 ) + V((a + 6)u t ) + e p (Qt - 0$) = f i - f 2 , (5.1) 
X t + BXi - SX 2 + /3(X!) - /3(X 2 ) + 7 (Xi) - 7 (X 2 ) (5.2) 

3 _ 6 ' (Xl) _ ^Hj _ _ 6-(x 2 )) £(U2)R 2 e£(U2) + e; - . 

Now, we test (|5.ip by u t and integrate in time. Recalling (|2.3p . it is not difficult to infer 

\\\Mt)\\ 2 L2 (n) +5 J llUtl^x^ds^ i||vJ-V 2 || 2 2(n) + j ||fl-f 2 ||H-i(n)||Uf|| ff i( a) ds + / 9 +i"lO+Ill, 

(5.3) 
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where we have 

h = - I (£(fe(Xi)u),u t ) Hl(n . Rd) ds<C J ||fo(Xi)|| L oo (f2) \\u\\ H i (n) ||u t || H1(n )ds 

S rt 



-lj !l u *llff 1 (r}) ds + C ' /„ ll u llHi(^) ds ' 



(5.4) 



whereas, the Lipschitz continuity of b on bounded intervals (cf. (|2. 18[) 1 and the Holder inequality yield 
ho = I (£((K X - &(X 2 ))u 2 ), u t ) H i (Qm ds 

J 

< C / ||u 2 || W / 1 .6(n)||X||L3 (f2 )|ju t || i /i ( n)d.s 

° rt (5-5) 



6 f f 
-iio H u *ll* 1 (n) dfl + C H Ua ll*-(o. r i^( n )) /„ l,X| 

(5 /"^ /"^ 
< 4 ^ Hu f ||?fi(n)ds + ^ J \\VX\\l* m &s + C J 



L 3 (o)ds 



IMIi 2 (n) ds 



where in the last inequality we have exploited the embeddings _ff 1 (J7) <<= L 3 (f2) C L 2 (il) and [521 Thm. 
16.4, p. 102], with v > a suitable constant to be chosen later and the constant C also depending on 

ll u 2lli- ( o,T;^ (n)) . Moreover, we get 

r* 

hi =P 



> f / (e* - e*)div( Ut )dxds 

Jo Jq 

<{f o Hu t |& 1( n)ds + P 2 C 5 ^ net - 9^||| 2(n) ds . 

Noting that ||u(i)||^ 1(0) 

dr ' we obtain from (|5.3[) (|5.5[) that 

1 S 

2\\Mt)\\h(Q) + I J l|ut||ffi W ds 

1 5 

< - v olli 2 (f2) + Cllfi - f 2ll^(0,T;ff-i(O)) + g y ll u dlHi(fi)ds + C\\ul - u 2 \\ 2 Hl{n) 

+ C f ( f S \\n t (r)\\ 2 H1{n) dr) ds + u [ || VX||| 2(n) d S + C f ||X|| 2 L2(0) ds 
Jo Vo / Jo Jo 

+ c P 2 f wei-e^l^ds. 

Jo 

Next, we test (|5.2|) by X integrate the resulting equation in time. With elementary computations, also 
taking into account the Lipschitz continuity of 7 (|2.20p . the monotonicity of j3 (|2.19[) . and the crucial 
inequality (|2.30[) . we get 

~IM*)ll£s(Q)+C9«^ ||VX||£ 2(fi) ds< \\\xl-Xl\\ 2 L 2 {nm +G ||X||| 2(fi) d S + / 12 +I 13 , (5.6) 
h 2 :=-!*[ 6'(X l} ( £(Ul)R f (Ul) - £(U2) y (U2) ) Xdxds, 



Jn 



'13 



((6'(xr) - ^(x 2 )) £(u2)R 2 e£(u2) x + (e; - es)x) dxd s . 
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Using (pTH|) and the fact that &'(Xi) G L°°(Q) we get 




where the last inequality is obtained arguing as for (|5.5I) . Again exploiting the Lipschitz continuity of b' 
on bounded intervals and the bound for 112 in L°°(0,T; Hq(CI; R d )), we get 




where the last estimate also follows from the continuous embedding 1 (SI) C L 6 (f2) and the Young 
inequality. Collecting now (|5.6p - (|5.8p . we arrive at 




The application of the standard Gronwall lemma gives immediately the desired continuous dependence 
estimate (|2.73p . □ 



Remark 5.1. If we replace the p-Laplacian (|2.25[) with the linear s-Laplacian (|2.34j) in the equa- 
tion for X, the continuous dependence estimate of Theorem [3] can be performed without assuming a 
to be constant (cf. (|2.72[0 . Indeed, in this case we would be able to deal with the additional term 
So (V((a(X 1 )-a(X 2 ))dfa 2 ,d t 

u )ff 1 (O K d ) d s > which results from subtracting the equations fulfilled by so- 
lution pairs (u;,Xj), i = 1,2. It would be possible to estimate it by means of the iJ s (fi)-norm of 
X = Xi — X2, which would pop in on the left-hand side of (|5.6[) . 

6 Proofs of Theorems |4] and [5] 
6.1 Proof of Theorem [4] 

Step 0: approximate equations. Equations (|3.17[) and (|3 . 18[) can be rephrased in terms of the interpolants 
w T , w T , u T , u T , u r , X T , X T , in a way analogous to (|4.ip - (|4.2[) . 
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Furthermore, taking into account that a = 0-?(_ao,o] an d that (3 = dI[ +oo j we observe that the 
minimum problem (|3.19[) yields for every k = 1 , . . . , K T that 



dx + $(r?) - 4>(X*) + / (7(r/)-7(X^) + ^- 1 (r / ^)) dx > 



for all ? ? e W-" 1 ^) with < 77 and 77 < X^ 1 a.c. in f2 , 

(recall the short-hand notation (|3.26p for h^T 1 ). Writing necessary optimality conditions for the above 
minimum problem, we infer 





77-X^ 1 


2 




1 


2/2 ( 


r 




T 





(d t X T (t)(r) - X T (t)) + d(x, VX T (i)) • V(r? - X T (i)) + 7(*t(*))0? - X T {t)) + h T (t)( v - X T {t)) 
for all t € [0, T] and all 77 € W 1,p (n) with < 77 and 77 < X T (t) a.c. in fi, 

where we have used the short-hand notation (cf. (|3.26[1 ) 

h T {t) := b (X T (i)J 0{w T {t)). 

Letting 77 = vip + X T (t) in (|6.1[) and dividing the resulting inequality by v > 0, we deduce that 
(d t X T (t)tp + d(x,WX T (t)) ■ Vip + -f(X T (t))ip + h T (t)ip^ dx > for all t G [0,T] 
and all <p € s.t. there exists v > with < i/p + X T {t) < X T (t) a.e. in O. 



dx > 
(6.1) 

(6.2) 
(6.3) 



Choosing ip = —d t X T (t) (observe that it complies with the constraint above, upon taking v = r), we 
therefore obtain 



(\d t X T {t)\ 2 + d{x, VX T (t)) • V(a t X T (i)) + 7 (X T (t))0 t X T (7i)) dx 



< 



(6.4) 



h T (t)d t X T (t)dx + / (7(X T (t))-7(X r (t)))ftX T (t)dar. 



Therefore, upon summing (|6.4[) over the index fc we deduce the discrete version of the energy inequality 
(|2T75)) for all < s < t < T, viz. 

ft At) 



t T (s) 



\d t X T \ A dxdr + ^{X T (t T (t)))+ / W(X T (t T (i)))da: 



< $(X r (t T ( S ))) + / W(X T (t T (s)))dx + 



L 2 (0,T;L 2 (O)) ) 



tr(*) 



t T (s) Jn 



d t X T -b'(X T ) 



e(u r )R e(u r ) 



9(w r ) I dxdr 



(6.5) 



+ CT 1 ' 2 II d t X T I 



where we have estimated the last term on the right-hand side of (|6.4p using that |j7(X r (i)) — 7(X t (<))||l2( t . l 
Ct 1 ' 2 11 a t x T 11 L 2 (o,T ; L 2 (f2)), thanks to the Lipschitz continuity of 7. 

S'iep 1 : compactness. In view of the a priori estimates from Proposition 13.101 we infer that there exist a 
vanishing subsequence (Tk)k and limit functions (w, u, X) such that convergences (|4.5[) . (|4.6|l . (|4.10|) - (|4.12p 
hold true as A: — > 00. Observe that (|4.6p in particular yields that X > and Xt < a.e. in £1 x (0,T). 
Arguing as in the proof of [271 Lemma 5.11], we now prove that 



X Tfc ^X in U'(fl,T;W 1 ' p (Sl)). 



(6.6) 



Indeed, [23 Lemma 5.2] gives a sequence ((p Tk )k C i p (0, T; W+' p (fi)) n L°°(Q) of test functions for 
fulfilling 



^ rfc ^ X in £ p (0, T; W 1,p (fi)), < <p Tk < X Tfc a.c. in Q x (0, T). 



(6.7) 
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Observe that the first of (|6.7I) and convergences (|4.6[) yield in particular 

X Tk ^X in L 2 (0,T;L 2 (O)). (6.8) 

We have 

c 7 [ [ |VX T t - VX| p dxds 



o Jn 



(6.9) 



< [ [ (d(x, VX T J - d(x, VX)) • V(X Tfc - X)dxds 
Jo Jn 

< / / d(x,WX Tk )-W(X Tk -cp Tk )dxds+ / / d(x, VX r J • V(^ Tfc - X)dxds 
Jo Jn Jo Jn 

- [ [ d(x, VX) • V(X Tfc - X)dxds = I 14 + hs + ho 
Jo Jn 

where the first inequality follows from (|2.26[) and the second one from elementary algebraic manipulations. 
Now, choosing tp := tp Tk — X Tk in f|6.3[) (which we are allowed to do thanks to (|6.7[1 ) and integrating in 
time, we obtain 

hi= / (d t X Tk + 7(X T J + h Tk ) (ip Tk -X Tk )dxds < C\\ip Tk -X Tk \\ L 2 {a T . L 2 {n)) -> as k -> oo, 
Jo Jn 

due to the bounds (pPT)) and (|X4U|) . and to (ffTTj) and ([6T5]) . We also have 

I-TibI < ||d(x, VX Tfc )|| iP ' (0iT;iP ' (f2)) ||V(^ rfc -X)||iP( ,T;iP(n)) < CllV^r, - X)\\lp(0,T;Lp(C1)) -> 0, 

where the second inequality follows from (|2.22l) and (|3.40l) . and the last passage is due to (|6.7[) . Taking 
into account that X r(c — X in L p (0,T; W 1,p (f2)) by (|4.9p . we also prove that 7x6 — > as fc — > oo. In this 
way, from (|6.9p we conclude (|6.6[) . Observe that, (|6.6[) combined with the bound (|3.40l) then yields (|4.9[) . 

5iep 2: passage to the limit. Arguing in the very same way as for the proof of Thm. [TJ it is possible to 
prove that (w, u, X) solve equations (|2.57|) and (|2.58[) . It now remains to prove the variational inequality 
(I2.76|) , together with (|2.77p . and the energy inequality (|2.78[) . As for the latter, it is sufficient to pass to 
the limit as k -> oo in (|S3)l . For this, we use convergences (|43]> . (g^) , ([4TT0]l - P~T2]l . ([4TT5|l . (|4TT5]I . as 
well as (|6.6[) . which in particular yields 

$(X r (t r (s))) -> $(X(s)) for a.a.s 6 (0,T). 

Clearly, the last term on the right-hand side of (|6.5[) tends to zero. Since the argument for (|2.76p ^ (|2.77l) 
is perfectly analogous to the one developed in the proof of [UJ Thm. 4.4], we refer the reader to [57] for 
all details and here just outline its main steps. Passing to the limit in (|6.3p as 4- with suitable test 
functions from [27j Lemma 5.2], we prove that for almost all t € (0, T) 

J {x t {t)Cp + d(ar, VX(t)) ■ V£ + 7 (*(t))0 + 6 / (X(t)) g(u(f))I ^ £(u( * )) y - 9(w(*))v) dx > 

for all ip G W.i' p (fi) with = 0} D {X(i) = 0}. 
From this, arguing as in the proof of |27l Thm. 4.4] we deduce that for almost all t G (0, T) 

J (x t (t)if + d(x, VX(t)) • + 7(X(*))V + ^(X(t)) £(u(t)) ^ e£(u(t)) y> - («;(*))¥>) dx 

> / f 7 (X(t)) + b , (X(t)) g(u(t))Re£(u(t)) -eKt))Vyda for all <p G VF^ p (r!). 

J{X(t)=o} V 2 / 
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Therefore, we take 

:= -3 {X=0] (x,t) ( 7 (*0M)) + b'(X(x,t)) £{uM) ^ £{uM) - e( W (x,t))Y 

for a.a. (x,t) G X (0,T), 

UrX_ i denoting the characteristic function of the set {X = 0}. From (|6.10p we deduce that, with this £ 
inequality (|2.76[) holds. Moreover, it is immediate to check that £ also complies with (|2.77l) . 

iSiep 3: strict positivity (|2.85[) of the temperature. Suppose that (|2.84p holds: the discrete strict positivity 
(13.24[) and convergences (|4.10p yield that, in the limit, w(x,t) > w (x) = _1 (j?o(ie)) for almost all 
(x,t) e!lx (0,T). Therefore, (|2l5)l ensues. □ 



6.2 Proof of Theorem [5] 

Step 1: compactness. For the interpolants u r , u T , u r , u r , X T , X r , £ T of the solutions (uji, Xj:, )^ 1 of 
the discrete Problcm l3.6l estimates ()3.36|1 — ()3.40f) and (|3.47p - (|3.50p hold. Therefore, standard strong and 
weak compactness results yield that there exist (u, X) fulfilling (|2.53l) - (|2.54[) and a subsequence | 
such that convergences (|4.5p and (|4.6I) hold. Moreover, estimates (|3.47[) ~ (|3.50[) also imply that X has the 
enhanced regularity (|2.89[) , and that 

X Tk ^*X in L°°(0,T; W 1+a ' p (n)) n W^°°(0,T; L 2 (fl)) for all < a < ±, 

X rfc , X Tfc , X Tfc ^X inZ^T;^'^)). 

Furthermore, there exist C and £ such that, possibly along a further subsequence, 

C Tfc -*C mL°°(0,T;L 2 (n)), (6.13) 

Pr k (X Tk )^*Z hiL°°(0,T;L 2 (fi)). (6.14) 

5tep 2: passage to the limit. Relying on convergences (|4.5[) . (|4.6p . and Q3.35ap . we take the limit of the 
discrete momentum equation (|3.18p . As for (|3.22p . we observe that, thanks to estimate (|3.47[) and the 
second of (|6.12p . there holds 

■B(X r J -> 'B(X) weakly* in L°°(0, T; L 2 (Q)) and strongly in L°°(0, T; W 1,P (Q)*). (6.15) 

Therefore, also taking into account f|4. 1 5[) we pass to the limit in (|3.22p and conclude (u, X, £, C) fulfill 
(|2.8ip . with Q(w) replaced by 0*. Furthermore, combining (|4.7p with (|6.14p we have 

limsup / / f3 Tk (X Tk )X Tk dxds < / / £Xdxds. 
k^oo Jq Jn Jo Jn 

Thanks to [TU1 Prop. 2.5, p. 27] we conclude that £ £ P{X) a.e. in f2 x (0,T). Finally, testing equation 
(I3.22p by dtX Tk and integrating in time, with calculations analogous to Q3.58P we find for all t <G [0,T] 

rtrW 



r^w r _ 
limsup / / d t X Tk C Tk dxds 
fc->oo Jo Jn 

/M*) f r 

< -liminf/ / \d t X Tk \ 2 dxds - liminf $(X Tfc (<)) -liminf / /3 Tk (X Tk {t)) 
k ->°° Jo Jn *-'" > °° fe ^°° Jn 

— liminf / / ^{X Tk )dtX Tk deeds — liminf / / h T d t X Tk dxds 
fe ->°° Jo Jn k ^°° Jo Jn k 

< - f I \X t \ 2 dxds - $(X(t)) - / P(X(t)) - / / 7 (X)X t dxds 
Jo Jn Jo Jo j n 

[ hdtXdxds = [ [ XtCdxds 
o Jn Jo Jn 
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(here h r is as in ([621), with ©*r in P la ce of Q(w T ), and h := b'(X) e{u)R 2 Mu) - 8*), where the last 
inequality is due to (|4.6[) and f)6. 12f) . the Mosco-convergence of ((3 Tk )r k to /?, the Lipschitz continuity of 
7, and (|4.15|) . The last identity follows from equation (|2.81|) . The aforementioned tokens of maximal 
monotone operator theory allow us to deduce from the fact that 

limsup / / d t X Tk ( Tk dxds < / / X t £dxds, 
fc->oo Jo Jn Jo Jq 

that ( £ a(X t ) a.e. in 17 x (0,T), which concludes the proof, n 

Remark 6.1. Indeed, in the proof of Theorem [5] the fact that a = <9/(_oo,o] has never been specifically 
used, therefore Thm. [5] extends to a maximal monotone operator a as in (|2.92[) . observing that, up to 
perturbing a with an affine function, it is not restrictive to suppose that 

e a(0) and S(0) = 0, whence a(x) > for all x £ K. 



7 Analysis of the degenerating system 

We now address the passage to the degenerate limit 8 4 in the full system (|2.57p - (|2.59[) . For technical 
reasons which will be clarified in Remark 17.51 later on, we focus on the irreversible case fi = 1, and neglect 
the thermal expansion term in the momentum equation, i.e. take p — 0. Furthermore, we confine the 
discussion to the case, in which, for 8 > the coefficients of both the elliptic operators in (|2.58[) are 
truncated, cf. Remark 17.21 below . In particular, we will take the functions a and b of the form 

a(X) = X, b(X) = X, and replace both coefficients by X + S. (7.1) 

Remark 7.1. The choice a(X) = 1 — X and b(X) = X in (|2.58[) and the truncation of both coefficients 
would lead to the momentum equation 

Utt + V((l-X + 6)u t ) + E((X + 5)u)=f inH- 1 (0;E d ), a.e. in(0,T) (7.2) 

for which the asymptotic analysis 8 4- would be less meaningful in the case of an irreversible evolution 
for X. For, starting from an initial datum Xo G W 1,p (fl) with max Ig j; Xq(x) < 1, we would have 
1 — X(x,t) > 1 — max^gjj Xo(x) > for all (x,t) £ fi x (0, T]. Hence, the limit 8 4- would not lead to 
elliptic degeneracy in (|7,2I) . 

Observe that the ensuing discussion can be suitably adjusted to the choice a(X) = X, b(X) = 1 — X, 
which is meaningful for phase transition models. 

Remark 7.2. It seems to us that both the coefficients a and b need to be truncated when taking the 
degenerate limit in the momentum equation. Indeed, on the one hand the truncation of a allows us to 
deal with the main part of the elliptic operator in (|2.58[) . On the other hand, in order to pass to the limit 
in the quadratic term on the right-hand side of (|2.59[) , we will also need to truncate b, cf. (|7.35l) later on. 

Theorem @] guarantees that for every 8 > there exists a triple (wg, us, Xs) as in (|2.52p - (l2.54l) fulfilling 
the enthalpy equation (|2.57|) with p = 0, the momentum equation 

d tt u 5 - div{{X + S)R v e{d t u s )) - div((X + <5)R e e(u 5 )) = f in i/- 1 (^;M d ), a.e. in (0,T), (7.3) 
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(where for later convenience we have dropped the operator notation (|2.5p ). as well as 



d t X s (x,t)<0 for a.a. (x,t) G CI x (0,T), (7.4) 
( + &(t) + 7(X«(t))) V + d(a:, X s (t)) ■ Vy>) dx < ^(MV^iM*)) + Q( Ws (t)fj 



for all ip G Wj' p (Q), for a.a. i G (0, T), 

(7.5) 



with &(*,*) = -J { x, =0} (M) ( 7 (X*(M)) + giH^lMH^)) _ (^(M))) + (7.6) 

for almost all (x, t) Sflx (0, T), (changing sign in (|2.76[) and recalling (|6.1ip ). and the energy inequality 
(12.78|) . As observed in Section [2T4l the family (ws, us, Xs)s then fulfills for all t G (0,T] the energy 
inequality 

w s (t)(dx) + 1 / |a t u 5 (t)| 2 da;+ / / |6» t X 5 | 2 dx dr + / a vis ((X 5 + 5)5 t u 5 , $iif) dr 



Ji2 



ioel((Xj(t) + 5)Ui(t),Ui(t))+$(X«(t))+ / W(X*(t))da; 



+ 2 / ' V °' 2 dX + k ael ^ X ° + ^ Uo ' + $(X °' ) + 



(7.7) 



f • dtiis dzdr + / / gdx. 
o Jo Jo 

First of all, following [35], in Prop. [731 below we deduce from equations (|2.57|) . (|7.3[) . and from (JT77J) 
some a priori estimates for the family (ws, us, Xs)s, independent of S > 0. 

Let us mention in advance that estimate (|7.1ip for (ws)s holds true only for the solutions (ws,us, Xs)s 
obtained through the time-discretization procedure of Section 13.31 Such solutions shall be referred to as 
approximable. Indeed, on the one hand, Remark 13.131 ensures that the discrete estimates p.39[) - (13.431) 
are valid with constants independent on 5: hence they are inherited by the approximable solutions (ws)s, 
yielding estimate (|7.11[) below. On the other hand, the calculations developed for the Fourth a priori 
estimate in Sec. 13.31 suggest that, in order to prove (|7.1ip for all weak solutions (ws)s to f|2 . 57[) . it would 
be necessary to test (|2.57p by ip = H(ws) with II as in (|3.68p . This is not an admissible choice due to 
the poor regularity of ws . Since we do not dispose of a uniqueness result for the irreversible full system, 
we cannot conclude (|7.1ip for all weak solutions (in the sense of Dcf. I2.13P (ws)s, and therefore we will 
restrict to approximable solutions. 

As it will be clear from the proof of Prop. l7T3l estimates (|7.8p - (|7.10p instead hold for all weak solutions 
(us,X s ). 

Proposition 7.3. Assume Hypotheses (I), (II), and (IV) with (3 = -Tm+oo); conditions (|2.47p - (|2.5ip on 
the data f, g,$o, u o ; v o ; Xo, and suppose that a, b are given by (|7.ip ■ Then, there exists a constant S > 
such that for all S > and for all (ws,us,Xs)s (approximable) weak solutions to the irreversible full 
system, the following estimates hold 

||w«5|jL^(0,T;L 1 (f2)) + ||c*tll,5||L°c(o,T;L 2 (n;Ill' 1 )) + ||Xa||£«( ,T;W 1 .*(fi))nff I (0,T;£ a (n)) C^) 

+ \\W(X 5 )\\ L ^ i0 . T . L i m <S, 
\\y/X + SR V e(dtUs)\\ L 2^ 0T . L 2^ a . R dxd- ) - ) + ||v / X _ +5R e e(u5)|| L oo(o i Tji 2 (n;K"*x'«)) ^ S> (7-9) 

\\dttUs\\L 2 (0,T;H- 1 (fl;R d )) < S, (7-10) 
\\ w s\\L r (0,T;W 1 ' r {n))nBV{[0,T};W 1 - r ' (n)») - S- (7-11) 
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Proof. Estimates (|7.8p - (|7.9p are straightforward consequences of the energy inequality (|7.7I) . taking into 
account that J n W(Xs(t)) dx > —C for a constant independent of t £ [0,T], estimating 

/ / f -9 t u 5 dxdr<i / / |f| 2 dx+i / / \d t u s \ 2 dxdr, 
Jo Jn z Jo Jn 1 Jo Jn 

and applying the Gronwall Lemma. Hence, (|7.10[) follows from a comparison in (|7.3p . in view of (|2.6a[) . 
Finally, (|7.1ip can be proved by observing that the discrete estimates (|3.42[i — (|3.43p are in fact independent 
of the parameter S > 0, hence they carry over to the approximable solutions (ws)s- □ 

As pointed out in [39] (see also [5]), estimates (|7.9D suggest that for the analysis S 1 it is meaningful 
to work with the quantities 



Us ■= VXs + Se(d t u s ), Vs ■= VXs + 6e{u s ) (7.12) 
in terms of which (|7.3p rewrites as 



d u us - div(^X s + SR V Mi ) - div(y/x s + 6R e r ls )=f in H'^Sl; R d ), a.e. in (0, T). (7.13) 

For later purposes, we also observe that, in the setting of (|7.ip and with notation (|7.12p . the total energy 
inequality (|2.79p for the triple (ws,us,Xs) can be reformulated as 



1 ' 



w s (t)(dx) + - I \d t u s (t)\ 2 dx + I I \d t X s \ 2 dxdr , 



\ [ T] s (t)R e T] s (t)dx + ^(Xs(t))+ j W(X s (t))dx 
^ Jn Jn 

< / w s (s)(dx) + l [ \d t u s (s)\ 2 dx + l [ r ls ( S )R e r ls ( S )dx + <£(Xs( S )) 



[ W(X s (s))dx+ [ [ f • d t u s dxdr+ [ [ gdxdr. 
Jn Js Jn J s Jn 



The following result shows that the limit <5 ^ preserves the structure (|7.13p of the momentum equation, 
as well as the enthalpy equation (|2.57p . The weak formulation (|2.76p - (|2.78p of the equation for X is 
generalized by ([7T"9" l) -([7^D ]l . cf. Rmk. El 

Theorem 6. Assume Hypotheses (I), (II), (IV) with (3 = -Zp.+oo)? an d (V) with <j> fulfilling (|2.26l) . 
Assume conditions (|2.47p - (|2.5ip on the data f, <?,$0i u o 5 v 0j -^Oj an d suppose that a, b are given by (|7.ip . 
Then, there exist w as in (|2.52p , and 

u e W hoo (0,T; L 2 (n;R d )) D H 2 (0,T; H- 1 ^;^)), /iG L 2 (0,T; L 2 (fi; M. dxd )), (7.15) 
r] e L°°(0, T; L 2 (tt] R dxd )), 

Xe L oo (0,T;PF 1 ' p (fi))niJ 1 (0,T; J L 2 (rj)), X(x,t) > 0, X t (x,t) < for a.a. (x, t) e ft x (0, T), 

(7.16) 

smc/i £/ia£ 

= \fXe(\i t ), rj = v^e(u) a.e. in any open set icSlx (0,T) s.t. X > a.e. m A, (7-17) 
fulfilling the weak enthalpy equation (I2.57P with p — 0, the weak momentum equation 

d tt u-div(\/XR v fj,)~dW{VxR e r]))={ in i? _1 (fi; R d ), a.e. in (0,T), (7.18) 
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as well as 

i-T r , pT 



[ [ ((d t X + 'y(X))ip + d(x,VX)- V(p^dxdt < ( ( (-^-r l R e r l + Q{w)\ipdxdt 
/or p e i p (0,T; W+ P (fi)) ni°°(Q) iwiift supp(^) C {X > 0}, 
and the total energy inequality for almost all t £ (0, T] 



(7.19) 



Jf(i) 



/ ( \X t \ 2 dxdr + \ ( ( /i(r)R„/i(r)da;dr 



< { w dx + \ [ |v | 2 dx+ia cl (X u ,u ) + $(X )+ / W(X )da: (7.20) 
Jn z Jn z Jn 

+ i ■ u t dxdr + J g dxdr 

Jo Jn Jo Jn 

withW(t)> [ w(t)(dx) + ^ [ \d t u(t)\ 2 dx + §(X(t)) + [ W(X(t))dx + d(t), 
Jn 2 J Jn 

where 3(t) := -liminf f r) Sk (t)R e r] Sk (t)dx, 
2 s k io J n 

(with (rjg ) a suitable subsequence of (t}§) from (|7.12[) ), and for all < t\ < t% < T there holds 
3i(r)dr> [ ( [ w(r)(dx) + $(X(r)) + / ( ^\d t u{r)\ 2 + W{X(r) + ^rj(r)R e r](r)\ dx ) 



(7.21) 



n Jn 



(7.22) 



Remark 7.4. Let us briefly compare the concept of weak solution (to the degenerating irreversible full 
system (|2.57[l - (|2.59[l ) arising from (j7.17p - (|7.20p , with the notion of weak solution (to the non- degenerating 
irreversible full system (|2.57p - (|2.59p ) given in Definition 12.131 in the case in which a(X) = b(X) = X. 
Suppose that the functions (u, X) in (|T. 15[) and (|7.16p have further regularity properties (|2.53[l - (|2.54|l . 
and that X > a.c. in O x (0,T). Then, (J7TT7J) holds a.e. in x (0,T), hence it is immediate to realize 
that (|7.19p coincides with (|2.76p . Furthermore, subtracting from (|7.20l) the weak enthalpy equation (|2.57l) 
tested by 1, we obtain a generalized form of the energy inequality (|2.78p for almost all t £ (0,T] and for 
s = 0. 

Proof. It follows from estimates (|7.8p - (|7.1ip and the same compactness arguments as in the proofs of 
Thins. [T] and [4] that there exist a vanishing sequence 8k I and functions w as in (|2.52l) and (u, X, H,r]) 
as in (|7.15p - (|7.16p such that as k — s- oo 

w Sk -> w in L v (0,T;W 1 - e ' r (Q.))nL a (0,T;L 1 (Q,)) for all e € (0,1] and all 1 < s < oo, (7.23) 

u Sk ^-*u in W hoo {0,T; L 2 (n;R d )) D H 2 (0,T; H-\n:R d )), (7.24) 

H Sk - n in L 2 (0, T; L 2 (il; M dxd )), (7.25) 

V Sk ^*V in L°°(0,T;L 2 (O;M dxd )), (7.26) 

X,5 fc ^*X in L oo (0,T;W 1 ' p (n))n J H" 1 (0,T;L 2 (Q)) ) (7.27) 

X 5fc ^XinC°([0,T];C o (n)), (7.28) 

the latter convergence due to the compactness results in [HI] and the compact embedding W 1,p (Sl) <s 
C°(U). Observe that (j?T2"4) and ([7T27|) respectively yield 

d t u 5k u t in C« cak ([0,T];i 2 (O;R d )), X Sk -+ X in C° eak ([0, T]; ^(Q)). (7.29) 
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From (|7.23|) . exploiting (|2.43[) we deduce that 

Q(w Sk )^e(w) in L 2 (0,T;L 2 (fi)). (7.30) 

Thus, we are in the position of passing to the limit as 5k i in (|2.57[) for the functions (ws k -,X$ k ), and 
conclude (|2~57|) for (w,X). 

Exploiting (|7.28[) and the fact that t n- X$(x,t) is nonincrcasing for all x € fi, with the very same 
argument as in the proof of [39l Prop. 4.3] it is possible to prove that fi and r\ have the form (|7.17p . In 
order to do that we can use the boundedness of e(u,5 fc ) and of e(d t us k ) in L 2 (K; R dxd ) for any compact 
cylinder K of the form Kq x [0, t] on which X > 0. Notice that on these cylinders X > 6 > for some <5 > 0. 
Hence, exploiting convergence (|7.28[) . we infer that there exists i5o > such that, for any < 6k < So, 
we have Xs k (x, t) + 5k > 5 for all x G K . Thus also Xs k (x, s) + 5k > 5 for all (x, s) e K = K x [0, t] 
because 1 i— > Xs(x, t) is nonincrcasing for all x € fi. Then we can identify at the limit fi and r\ and cover 
A in (|7.17j) by cylinders of the form K above. Hence, relying on (|7.24p - (|7.28j) it is immediate to pass to 
the limit in (|7.13[) and conclude (|7.18l) . 

Next, we prove that 

X 3k ^X mL p (0,T;W hp (n)). (7.31) 

For this, we repeat the arguments from Step 2 in the proof of Thm. 0] based on [37J Lemma 5.11]. 
Namely, we apply [27J Lemma 5.2], which gives a sequence (tps k )k C L p (0,T;W^ p (fl)) (1 L°°(Q), such 
that tp Sh -> X in 2?(0,T; W 1,p (fi)) and < tp Sh < X Sk a.e. in fi x (0,T). Relying on assumption (|2~5d) . 
with the same calculations as in (|6.9p we then have 

c 7 / / \VX Sk - VX| p dxds 

< f [ d{x,VX 5k )-V{X Sk -<p 5k )dxds + [ [ d{x,VX Sk )-V{<p Sh -X)dxds (7.32) 
Jo Jn Jo Jn 

-I f d(x,VX)-V(X Sk -X)dxds = I 17 + I 18 + I 19 . 
Jo Jn 

Now, choosing ip$ k := ip$ k — Xg k as a test function for (|7.5[) and integrating in time, we obtain 

/it = f [ (d t X Sk + £ Sk + j(X Sk ) + £ ( u ^)y( u ^) _ Q{wSk) ) {lfSk - XSk )dxds = J? 7 + I\ 7 + Jf 7 , 
where 



o Jn 



I vi = I I ( d t*s k +l{X Sk ) - &(w Sk )) (ip Sk -X Sk )dxds < C\\(p Sk - X 5 J| L 2 (0jT . i2(n)) -> as fc -> oo, 



J i7 = / (&*> a-^JffMfn) ds < 

f £KJ y (uA " J (^-x A -jdxd.s< 



the second inequality due to (|2.77p . and the third one to the fact that ips k < X,5 fc a.e. in fi x (0,T). 
Calculations completely analogous to the ones developed in the proof of Thm. U yield that Jig — s- 
as k — > oo. In this way, from (|7.32p we conclude (|7.31l) . 

We are now in the position to pass to the limit in (the time-integrated version of) (|7.5p and conclude 
(|7.19p . To this aim, we observe that, for any fixed test function <p as in (|7.19l) . supp((p) is a compact 
subset of fi x [0,T]. Hence there exists X > such that X(x,t) > X > for all (x,t) <E supp(y), and, by 
(|7.28p . there exists k e N such that for k > k 

Xs k (x, t)>^X>0 for all (x, t) e supp(y). (7.33) 
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Therefore, 



k- 



o Jn 



lim / / fatpdxdt = 0, (7.34) 



since supp(£,5) C {Xg = 0} by ()7.6[l . Also exploiting (|7.31[) . we succeed in taking the limit of the left-hand 
side of (|7.5|) . As for the right-hand side, we use (|7.30|) and argue in the following way 

e(u Sk )R e e(u Sk ) , \ _ ,. . . f T f 1 




2 



-^dzdt ] = -hminf / / — ——-T] g R e r] g ipdxdt 

I s k ^o J Q J n 2(X Sk +S k ) 

' t (7.35) 

< — / / — n R e riipdxdt , 

- y y n 2x 

where we have used that, thanks to (|7.28|) and (|7.33|) . 2 (X^ +s .) ~^ IX un if° rm ly on su PP(v)i thus the 
last inequality e.g. follows from the lower semicontinuity result of [2]. 

Finally, (|7.20[) follows from taking the limit as 5k — > of the total energy inequality (|7.14[) , written on 
the interval (0,t) for any t € (0, T]. Observe that, by (|7.14p (cf. also the arguments in the proofs of [3"51 
Prop. 4.3]), the map 

t^K s (t) := [ w s {t)dx + l [ \d t u s (t)\ 2 dx + \ [ r) s (t)R e ri s (t)dx + $(Xs(t))+ [ W(X s (t))dx 
Jn 1 Jn 2 Jq Jq 

has (uniformly) bounded variation. Therefore, by Helly's theorem up to a subsequence there exists !K 
such that !Ks k (t) —> "K{t) for all t £ [0, T]. To identify J{, we take the liminf as Sk —> of the first, 
second, fourth, and fifth term in CK<5 (t), exploiting convergences (|7.23[) . (|7.27p . (|7.28p . (|7.29|) . as well as 
(|7.3ip . and relying on lower semicontinuity arguments. Therefore we conclude that (|7.21j) holds. Finally, 
inequality (|7.22[) follows combining the following facts: on the one hand, since ('Ks)g C L°°(Q, T) is 
uniformly bounded due to estimates (|7.8p - (|7.1ip . the dominated convergence theorem ensures 

/■*2 rt2 

/ , K Sk (r) dr -> / JC(r) dr as k -> oo for all < h < t 2 < T. 

Jt x Jt x 

On the other hand, on account convergences (|7.23p ~ (|7.28p . by weak lower semicontinuity arguments we 
have that lim inf fc^oo J t * 2 Jis k (r) dr is greater or equal than the right-hand side of (|7.22p . This concludes 
the proof. □ 

Remark 7.5. As it is clear from the above lines, the proof of Thm. [5] strongly relies on the following 
properties: 

1. the compact embedding of VF 1,P (57) into C°(f2); 

2. the fact that t \— > X$(t,x) is nonincreasing for all x £ f2, which follows from the irreversibility 
constraint. 

These are the reasons why we have restricted the analysis of the degenerate limit to the irreversible 
system. Within this setting, we further need to assume p = 0. Indeed, because of the lack of estimates 
on div(u t ) for 6 J, 0, we would not be able to the limit in the term pdiv(u t )<d(w) in (|2.57p as 6 I 0. 

We also point out that, seemingly, the total energy inequality (|7.20l) cannot be improved to an inequality 
holding on any subinterval (s, t) C (0, T). Indeed, for the sequence (r)g )k only the weak convergence 
(|7.26p is available, which docs not allow us to take the limit of the right-hand side of (|7.14p but for s = 0. 

Finally, observe that the proof of Thm.[6]simplifies if the operator 25 is given by the nonlocal s-Laplacian 
operator A s . In this case, in order to pass to the limit in (|7.5|l it is no longer necessary to prove the 
strong convergence (|7.3ip for (X$ k )k- In fact, the term d(x, X$ k ) ■ V(p in (|7.5p is replaced by a s (Xs k7 f), 
which can be dealt with by weak convergence arguments due to the linearity of the operator A s . 
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